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ABSTRACT 

Aims.  We  present  the  first  multi-epoch  study  that  includes  concurrent  mid-infrared  and  radio  interferometry  of  an  oxygen-rich  Mira 
star. 

Methods.  We  obtained  mid-infrared  interferometry  of  S  Ori  with  VLTI/MIDI  at  four  epochs  in  December  2004,  February/March 
2005,  November  2005,  and  December  2005.  We  concurrently  observed  v  =  1,  J  =  1-0  (43.1  GHz),  and  v  =  2,  J  =  1-0  (42.8GHz) 
SiO  maser  emission  toward  S  Ori  with  the  VLBA  in  January,  February,  and  November  2005.  The  MIDI  data  are  analyzed  using 
self-excited  dynamic  model  atmospheres  including  molecular  layers,  complemented  by  a  radiative  transfer  model  of  the  circumstellar 
dust  shell.  The  VLBA  data  are  reduced  to  the  spatial  structure  and  kinematics  of  the  maser  spots. 

Results.  The  modeling  of  our  MIDI  data  results  in  phase-dependent  continuum  photospheric  angular  diameters  of  9.0  ±0.3  mas 
(phase  0.42),  7.9±0.1mas  (0.55),  9.7±0.1mas  (1.16),  and  9.5±0.4mas  (1.27).  The  dust  shell  can  best  be  modeled  with  A1203 
grains  using  phase-dependent  inner  boundary  radii  between  1.8  and  2.4  photospheric  radii.  The  dust  shell  appears  to  be  more  compact 
with  greater  optical  depth  near  visual  minimum  (tv  ~  2.5),  and  more  extended  with  lower  optical  depth  after  visual  maximum 
( Ty  ~  1.5).  The  ratios  of  the  43.1  GHz/42.8  GHz  SiO  maser  ring  radii  to  the  photospheric  radii  are  2.2  ±  0.3/2. 1  ±0.2  (phase  0.44), 
2.4  ±  0.3/2. 3  ±  0.4  (0.55),  and  2.1  ±  0.3/1. 9  ±  0.2  (1.15).  The  maser  spots  mark  the  region  of  the  molecular  atmospheric  layers  just 
beyond  the  steepest  decrease  in  the  mid-infrared  model  intensity  profile.  Their  velocity  structure  indicates  a  radial  gas  expansion. 
Conclusions.  S  Ori  shows  significant  phase-dependences  of  photospheric  radii  and  dust  shell  parameters.  ALO3  dust  grains  and 
SiO  maser  spots  form  at  relatively  small  radii  of  ~1.8-2.4  photospheric  radii.  Our  results  suggest  increased  mass  loss  and  dust 
formation  close  to  the  surface  near  the  minimum  visual  phase,  when  ALO3  dust  grains  are  co-located  with  the  molecular  gas  and  the 
SiO  maser  shells,  and  a  more  expanded  dust  shell  after  visual  maximum.  Silicon  does  not  appear  to  be  bound  in  dust,  as  our  data  show 
no  sign  of  silicate  grains. 

Key  words,  techniques:  interferometric  -  masers  -  stars:  AGB  and  post-AGB  -  stars:  atmospheres  -  stars:  mass-loss  - 
stars:  individual:  S  Orionis 


1.  Introduction 

The  evolution  of  cool  luminous  stars,  including  Mira  variables, 
is  accompanied  by  significant  mass  loss  to  the  circumstellar  en¬ 
vironment  (CSE)  with  mass-loss  rates  of  up  to  1 0  4  M0/ year 


*  Based  on  observations  made  with  the  Very  Large  Telescope 
Interferometer  (VLTI)  at  the  Paranal  Observatory  under  program 
IDs  074.D-0075  and  076.D-0338. 

**  Based  on  observations  made  with  the  Very  Long  Baseline  Array 
(VLBA)  under  project  BB192.  The  VLBA  is  operated  by  the 
National  Radio  Astronomy  Observatory  (NRAO).  The  National  Radio 
Astronomy  Observatory  is  a  facility  of  the  National  Science  Foundation 
operated  under  cooperative  agreement  by  Associate  Universities,  Inc. 

***  Tables  A.  1-A.4  are  only  available  in  electronic  form  at  the  CDS  via 
anonymous  ftp  to  cdsarc.u-strasbg.fr  (130.79.128.5)  or  via 
http: //cdsweb.u-strasbg . fr/cgi-bin/qcat?J/A±A/470/191 
t  Color  versions  of  Figs.  2-8,  10,  and  11  are  available  in  electronic 
form  via  http :  //www.  aanda.  org 


(e.g.,  Jura  &  Kleinmann  1990).  This  mass-loss  process  signif¬ 
icantly  affects  any  further  stellar  evolution  and  is  one  of  the 
most  important  sources  for  chemical  enrichment  of  the  interstel¬ 
lar  medium.  The  detailed  nature  of  the  mass-loss  process  from 
evolved  stars,  and  especially  its  connection  with  the  pulsation 
mechanism  in  the  case  of  Mira  variable  stars,  is  a  matter  of  cur¬ 
rent  investigation. 

The  conditions  close  to  the  stellar  surface  can  be  studied  well 
by  means  of  optical  long-baseline  interferometry.  This  technique 
has  provided  information  regarding  the  stellar  photospheric  di¬ 
ameter,  effective  temperature,  center-to-limb  intensity  variation 
(CLV),  and  atmospheric  molecular  layers  for  a  number  of  Mira 
variables  (see,  e.g.,  Quirrenbach  et  al.  1992;  Haniff  et  al.  1995; 
van  Belle  et  al.  1996;  Thompson  et  al.  2002;  Mennesson  et  al. 
2002;  Hofmann  et  al.  2002;  van  Belle  et  al.  2002;  Perrin  et  al. 
2004;  Millan-Gabet  et  al.  2005). 

Woodruff  et  al.  (2004)  and  Fedele  et  al.  (2005)  have  re¬ 
cently  shown  that  observed  near-infrared  Af-band  visibilities  of 
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the  oxygen-rich  prototype  Mira  variables  o  Cet  and  R  Leo  are 
very  different  from  uniform  disc  (UD)  models  already  in  the  first 
lobe  of  the  visibility  function,  and  that  they  correspond  closely 
to  predictions  by  self-excited  dynamic  Mira  model  atmospheres 
that  include  effects  from  molecular  layers  (models  by  Hofmann 
et  al.  1998;  Tej  et  al.  2003b;  Ireland  et  al.  2004a,b).  Ohnaka 
et  al.  (2006b)  studied  the  comparison  of  the  same  dynamic  Mira 
model  atmospheres  with  mid-infrared  interferometric  and  spec¬ 
troscopic  observations.  Recently,  Ireland  &  Scholz  (2006)  added 
the  formation  of  dust  in  a  self-consistent  way  to  the  same  dy¬ 
namic  Mira  model  atmospheres.  They  find  that  dust  would  form 
at  2-3  times  the  average  photospheric  radius  for  certain  plausible 
parameter  values. 

The  structure  of  the  atmospheric  molecular  shells  located 
above  the  photosphere,  as  well  as  the  dust  shell,  can  be  probed  by 
mid-infrared  interferometry.  This  has  been  successfully  demon¬ 
strated  by  Ohnaka  et  al.  (2005)  using  the  spectro-interferometric 
capabilities  of  the  VLTI/MIDI  facility  for  observations  of  the 
oxygen-rich  Mira  star  RR  Sco.  The  model  used  in  this  work 
includes  a  warm  molecular  layer  consisting  of  SiO  and  H2O, 
as  well  as  an  optic  ally-thin  dust  shell  of  ALO3  and  silicate. 
Recent  VLTI/MIDI  observations  of  the  carbon-rich  Mira  star 
V  Oph  by  Ohnaka  et  al.  (2007)  indicate  that  carbon-rich  Miras 
also  have  extended  atmospheric  layers  of  polyatomic  molecules 
(C2H2)  and  dust  shells  (amorphous  carbon  and  SiC).  Information 
on  dust  shells  around  Mira  variables  has  also  been  derived 
using  mid-infrared  interferometry  with  the  Berkeley  Infrared 
Spatial  Interferometer  ISI  (see  e.g.,  Danchi  et  al.  1994;  Weiner 
et  al.  2006;  Tatebe  et  al.  2006).  Ireland  et  al.  (2005)  recently 
used  0.9 /rm  interferometric  polarimetry  measurements  of  the 
Mira  variables  R  Car  and  RR  Sco  to  place  constraints  on  the 
distribution  of  the  dust  shell. 

Complementary  information  regarding  the  molecular  shells 
around  oxygen-rich  AGB  stars  can  be  obtained  by  observing  the 
maser  radiation  that  some  of  these  molecules  emit.  Maser  emis¬ 
sion  from  the  three  most  common  maser  molecules,  SiO,  H2O, 
and  OH,  traces  regions  of  the  CSE  on  scales  from  a  few  to  sev¬ 
eral  hundred  AU.  Masers  provide  a  unique  probe  of  structure, 
kinematics,  and  polarization  properties  of  the  environment  of 
these  stars.  The  SiO  maser  radiation  in  the  CSE  of  oxygen-rich 
Mira  variables  has  been  mapped  using  the  Very  Long  Baseline 
Array  (VLBA)  by,  e.g..  Diamond  et  al.  (1994),  Kemball  & 
Diamond  (1997),  Boboltz  et  al.  (1997),  Diamond  &  Kemball 
(2003),  Cotton  et  al.  (2004),  Boboltz  &  Wittkowski  (2005,  here¬ 
after  BW05),  Boboltz  &  Diamond  (2005),  and  Cotton  et  al. 
(2006).  The  SiO  maser  emission  is  typically  found  to  arise  from 
a  clumpy  ring  within  a  few  stellar  radii  of  the  photosphere,  in¬ 
dicating  a  tangential  amplification  process.  Humphreys  et  al. 
( 1996,  2002)  have  predicted  the  location  at  which  SiO  molecules 
form  in  the  CSE  of  Mira  variables  by  combining  maser  models 
with  hydrodynamic  pulsation  models.  Their  model  predictions 
are  roughly  consistent  with  the  measurements  mentioned  above. 

Observational  results  regarding  the  detailed  relationships  be¬ 
tween  the  stellar  photosphere,  the  molecular  layer,  the  dust  shell, 
and  the  SiO  maser  ring  often  suffer  from  uncertainties  inherent 
in  comparing  observations  of  variable  stars  widely  separated  in 
time  and  stellar  phase,  as  discussed  in  BW05.  To  overcome  this 
limitation,  we  have  established  a  program  of  concurrent  infrared 
interferometry  using  the  VLTI  and  radio  interferometry  at  the 
VLBA.  The  former  aims  at  constraining  the  photospheric  radius, 
the  characteristics  of  atmospheric  molecular  layers,  and  parame¬ 
ters  of  the  dust  shell.  The  latter  aims  at  concurrently  mapping  the 
SiO  maser  emission.  Our  final  goal  is  a  better  understanding  of 
the  mass-loss  process  and  its  connection  to  stellar  pulsation.  Our 


pilot  study  in  BW05  on  the  Mira  variable  S  Ori  included  coordi¬ 
nated  near-infrared  A'-hand  interferometry  to  constrain  the  stel¬ 
lar  photospheric  diameter  and  VLBA  mapping  of  the  SiO  maser 
radiation  toward  this  source. 

A  further  uncertainty  in  comparing  photospheric  radii  to  the 
extensions  of  the  dust  and  maser  shells  arises  from  the  compli¬ 
cation  that  near-  and  mid-infrared  CLVs  of  finite  bandwidth  in¬ 
clude  a  blend  of  intensities  from  continuum-forming  layers  and 
overlying  molecular  layers.  This  effect  has  often  resulted  in  over¬ 
estimated  continuum  photospheric  diameter  values  (see,  e.g.  the 
discussions  in  Jacob  &  Scholz  2002;  Mennesson  et  al.  2002; 
van  Belle  et  al.  2002;  Woodruff  et  al.  2004;  Ireland  et  al.  2004a,b; 
Perrin  et  al.  2004;  Fedele  et  al.  2005).  An  overestimated  photo¬ 
spheric  diameter  would  result  in  biased  relative  distances  of  the 
dust  shell  and  the  maser  ring  from  the  stellar  surface,  even  if  ob¬ 
tained  at  the  same  stellar  phase  and  cycle.  A  detailed  comparison 
of  observations  to  dynamic  model  atmospheres,  as  mentioned 
above,  can  be  used  to  relate  the  observable  quantities  to  the  con¬ 
tinuum  photospheric  radius  and  thus  to  overcome  this  limitation. 

Here,  we  present  VLTI/MIDI  mid-infrared  interferometry 
of  S  Ori  at  4  epochs/stellar  phases  and  coordinated  VLBA 
mapping  of  the  SiO  42.8  GHz  and  43.1GHz  maser  transitions 
at  3  epochs/stellar  phases  that  are  contemporaneous  to  the  first  3 
of  our  4  MIDI  epochs.  The  first  two  epochs  are  located  near  the 
stellar  minimum,  and  the  later  epochs  shortly  after  the  following 
stellar  maximum. 


2.  Lightcurve  and  characteristics  of  S  Ori 

S  Ori  is  a  Mira  variable  star  with  spectral  type  M6.5e-M9.5e 
and  V  magnitude  7.2-14.0  (Samus  et  al.  2004).  Templeton  et  al. 
(2005)  report  that  S  Ori’s  period  appears  to  vary  in  a  seem¬ 
ingly  sinusoidal  fashion  between  about  400  and  450  days  over 
about  the  past  100  years.  The  mean  period  is  P  =  420.0  days 
with  a  latest  period,  corresponding  to  our  observation  epochs,  of 
P  ~  430  days.  We  find  a  good  correspondence  with  the  AAVSO 
(Henden  et  al.  2006)  and  AFOEV  (CDS)  data  for  the  cycle  of 
our  observations,  as  well  as  for  the  most  recent  ~5  cycles  using 
P  —  430  days  and  a  Julian  Date  of  last  maximum  brightness 
To  =  2453  190  days.  Figure  1  shows  the  visual  lightcurve  of 
S  Ori  as  a  function  of  Julian  Date  and  stellar  phase  based  on 
these  values.  Also  indicated  are  the  dates/epochs  of  our  VLTI 
and  VLBA  observations  as  described  in  Sects.  3  and  4. 

The  distance  to  S  Ori,  as  for  many  Mira  stars,  is  not  well 
known.  Van  Belle  et  al.  (2002)  have  estimated  the  distance  to 
S  Ori  to  480  pc  +  120  pc  based  on  a  calibration  of  the  period- 
luminosity  relationship  by  Feast  et  al.  (1989),  and  we  adopt  this 
value. 

The  near-infrared  K-band  UD  angular  diameter  ©fiD  of  S  Ori 
has  been  measured  by  van  Belle  et  al.  (1996),  Millan-Gabet  et  al. 
(2005),  and  by  BW05  to  values  between  9.6  mas  and  10.5  mas  at 
different  phases.  Diameter  measurements  of  S  Ori  at  other  wave¬ 
lengths  have  so  far  not  been  reported. 

S  Ori  exhibits  SiO  and  OH  maser  emission  (Benson  et  al. 
1990),  while  a  detection  of  H2O  maser  toward  S  Ori  has  not 
been  reported.  The  maps  of  the  SiO  maser  emission  in  BW05 
are  the  first  toward  S  Ori.  They  derived  average  distances  and 
standard  deviations  of  the  SiO  maser  spots  from  the  center  of 
their  distribution  at  phase  0.73  of  9.4+  1.4  mas  and  8.8  +  1.7  mas 
for  the  43.1  GHz  and  42.8  GHz  transitions,  respectively. 

Sloan  &  Price  (1998)  report  on  a  relatively  low  dust- 
emission  coefficient  (Dec,  the  total  emission  of  the  dust  to  the 
total  emission  of  the  star  in  the  wavelength  range  7.7-14.0  /m 1) 
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Table  1.  VLTI/MIDI  observations  of  S  Ori. 


Ep. 

Date 

Time 

UTC 

JD 

-2.4E6 

®viS 

Conf. 

B 

Disp. 

Elem. 

BC 

Bv 

[m] 

PA 

[deg] 

Seeing 

["] 

To 

[msec] 

A 

2004-12-31 

02:43 

53370.6 

0.42 

UT3-UT4 

62  m 

Prism 

HS 

62.27 

-71.9 

0.76 

3.4 

B 

2005-02-21 

02:37 

53422.6 

0.54 

UT2-UT4 

89  m 

Prism 

HS 

74.31 

83.0 

0.74 

3.6 

B 

2005-03-01 

02:39 

53430.6 

0.56 

UT3-UT4 

62  m 

Prism 

HS 

41.65 

-55.6 

0.62 

3.6 

B 

2005-03-03 

00:59 

53432.5 

0.56 

UT2-UT4 

89  m 

Prism 

HS 

84.35 

82.7 

0.98 

3.7 

C 

2005-11-12 

08:57 

53686.9 

1.16 

UT1-UT4 

130m 

Prism 

SP 

123.02 

62.8 

0.54 

4.4 

C 

2005-11-15 

07:21 

53689.8 

1.16 

UT2-UT3 

Aim 

Prism 

SP 

46.26 

43.9 

0.61 

2.9 

D 

2005-12-30 

03:00 

53734.6 

1.27 

UT1-UT4 

130m 

Prism 

SP 

123.89 

59.4 

1.12 

2.8 

a.  -5 


O) 

E 


b. 


-4 


-3 


Stellar  cycle/phase 


-2 


-1 


Julian  Date  -  2450000 


0.0 


Stellar  phase 


0.5 


1.0 


1.5 


Julian  Date  -  2450000 


Fig.  1.  Visual  lightcurve  of  S  Ori  as  a  function  of  Julian  Date  and  stel¬ 
lar  cycle/phase.  Data  are  from  the  AAVSO  (Henden  et  al.  2006)  and 
AFOEV  (CDS)  databases.  We  adopt  a  period  of  P  =  430  days  and  JD  of 
last  maximum  brightness  Tq  =  2453190.  The  arrows  indicate  the  dates 
of  our  VLTI  and  VLBA  observations,  a)  The  most  recent  5-6  cycles  il¬ 
lustrating  the  correspondence  of  the  data  with  our  adopted  P  and  To  val¬ 
ues.  The  indicated  dates  of  observations  are  the  epochs  from  BW05 
(JD  between  2600  and  2800)  and  the  epochs  of  the  present  paper  (JD 
between  3300  and  3800).  b)  Enlarged  view  covering  our  VLTI/MIDI 
and  VLBA  observations  presented  in  this  paper.  The  dates  of  observa¬ 
tions  are  combined  into  4  epochs  A,  B,  C,  D. 


for  S  Ori  of  0.24,  i.e.  similar  to  RR  Sco  (Dec  =  0.21)  or  R  Leo 
(Dec  =  0.23),  indicating  an  optically  thin  dust  shell. 


3.  VLTI/MIDI  measurements 

We  obtained  mid-infrared  (8-1 3 gm)  interferometry  of  S  Ori 
with  the  instrument  MIDI  (Leinert  et  al.  2003)  at  the  ESO  VLT 


Interferometer  (VLTI;  Glindemann  et  al.  2003)  in  service  mode 
between  31  December  2004  and  30  December  2005.  Light  was 
combined  from  two  8  m  Unit  Telescopes  (UTs).  We  chose  to  use 
the  prism  with  spectral  resolution  R  =  A/AA  =  30  to  disperse  the 
interferograms.  Two  various  techniques  of  beam  combinations 
were  used,  namely  HIGH_SENS  (HS,  2004-12-31  to  2005-03-03) 
and  SCI_PH0T  (SP,  2005-11-12  to  2005-12-30)  combinations. 
Using  the  HS  combination,  all  arriving  light  is  first  combined  in- 
terferometrically,  and  the  photometric  spectrum  for  each  incom¬ 
ing  beam  is  recorded  sequentially.  Using  the  SP  combination, 
beam  splitters  are  used  to  record  the  interferograms  and  the  pho¬ 
tometric  spectra  simultaneously.  An  estimate  of  the  photometric 
count  rates  at  the  time  of  recording  the  interferometric  data  is 
needed  to  compute  the  raw  visibility  values,  and  the  SP  tech¬ 
nique  is  thus  expected  to  result  in  a  higher  precision  of  the  visi¬ 
bility  measurement  for  bright  targets.  Our  observations  followed 
the  standard  procedures  for  MIDI  service  mode  observations  as 
described  in  detail  in  the  MIDI  user  manual1.  For  principles  of 
observations  with  MIDI  see  also  Leinert  et  al.  (2004). 

The  details  of  our  observations  are  listed  in  Table  1 .  Listed 
are  the  epoch,  date,  time,  and  Julian  Date  (JD)  of  the  obser¬ 
vation;  the  visual  phase  4>vis;  the  baseline  configuration  used 
with  its  ground  length  //;  the  dispersive  element  and  beam  com¬ 
biner  BC  (HS  for  High_Sens  and  SP  for  Sci_Phot);  and  the 
projected  baseline  length  Bp  and  its  position  angle  PA  on  sky 
(deg.  east  of  north).  The  last  two  columns  describe  the  ambi¬ 
ent  conditions  DIMM  seeing,  and  coherence  time  to,  both  at 
500  nm.  Dates  of  observation  were  combined  into  4  epochs  A, 
B,  C,  D.  Each  epoch  has  a  width  of  less  than  17  days,  i.e.  less 
than  5%  of  S  Ori’s  variability  period.  The  dates  of  observation 
and  the  epochs  compared  to  the  visual  lightcurve  are  indicated 
in  Fig.  1 .  We  estimate  the  uncertainty  in  assigning  a  visual  phase 
value  to  each  of  our  observing  epochs  to  ~0.1,  taking  into  ac¬ 
count  both  the  time  widths  of  our  epochs  and  the  uncertainty 
in  the  choice  of  S  Ori’s  current  variability  period  and  date  of 
last  maximum  (see  Sect.  2).  Epochs  A  and  B  are  located  around 
stellar  minimum  with  a  separation  of  about  2  months  (phase  dif¬ 
ference  ~0. 1 ),  while  epochs  C  and  D  occurred  about  a  year  later 
shortly  after  the  following  stellar  maximum  with  a  separation 
of  6  weeks.  Additional  MIDI/UT  data  on  S  Ori  from  2004-12- 
29,  2005-01-02,  2005-02-28,  2006-03-12,  2006-03-12,  as  well 
as  MIDI/ AT  data  taken  between  October  2005  and  March  2006 
were  not  used  because  of  various  technical  problems  with  these 
data  sets.  The  interferometric  field  of  view  of  MIDI  used  with 
the  UTs  is  ~250  mas,  which  corresponds  to  about  25  times  the 
expected  photospheric  angular  diameter  of  S  Ori. 

We  used  6  CMa  (HD  50778)  as  the  interferometric  calibra¬ 
tor,  a  K4  giant  at  a  distance  on  sky  of  27°,  which  was  observed 
within  about  30  min  of  each  S  Ori  observation.  A  number  of 


1  http : //www. eso . org/instruments/midi 
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Table  2.  Properties  of  the  MIDI  calibration  stars  used. 


HD 

Name 

Sp.T. 

Old 

[mas] 

Teff 

[K] 

F 12 

[Jy] 

( 

[°] 

50778 

6  CMa 

K4III 

3.90  ±  0.02 

4049 

24.6 

27 

31421 

o~  Ori 

K2III 

2.65  ±  0.01 

4500 

13.6 

12 

49161 

17  Mon 

K4III 

2.44  ±  0.01 

4248 

7.2 

20 

61935 

ff  Mon 

G9III 

2.24  ±  0.01 

4813 

10.2 

36 

94510 

u  Car 

K1III 

2.23  ±  0.01 

4951 

6.9 

90 

107446 

eCru 

K3.5  III 

4.43  ±  0.02 

4139 

32.4 

101 

observations  of  other  MIDI  calibration  stars  using  the  same  in¬ 
strument  settings  were  obtained  during  some  of  our  observation 
nights  for  other  programs,  and  we  make  use  of  these  data  as  well. 
Table  2  lists  the  calibration  stars  and  their  properties.  These  pa¬ 
rameters  are  taken  from  the  instrument  consortium’s  catalog2. 
Listed  are  the  HD  number,  the  name,  the  spectral  type,  the  limb- 
darkened  diameter  ©ld,  the  effective  temperature  Tes,  the  flux 
at  1 2  //m,  and  the  distance  on  sky  from  S  Ori  /.  The  UD  diam¬ 
eter  for  K  giants  in  the  mid-infrared  based  on  spherical  model 
atmospheres  is  smaller  by  ~5%  than  ©ld  (estimated  with  proce¬ 
dures  as  used  in  Wittkowski  et  al.  2006).  HD  50778  is  our  main 
calibrator.  The  other  calibrators  were  observed  during  some  of 
our  observation  nights  for  other  programs.  Among  the  interfer¬ 
ometric  calibration  stars  listed  in  Table  2,  absolutely  calibrated 
spectra  are  available  for  HD  31421  and  HD  61935  (Cohen  et  al. 
1999). 

3.1.  MIDI  data  reduction 

Mean  coherence  factors  and  photometric  count  rates  were  ob¬ 
tained  for  each  data  set  of  S  Ori  and  the  calibration  stars  using  the 
MIA+EWS  software  package,  version  1.4  (Jaffe,  Koehler,  Cotton, 
Chesneau,  et  al.3),  complemented  by  the  tool  MyMidiGui,  ver¬ 
sion  1.41  (Hummel4).  Data  for  32  spectral  channels  within  the 
nominal  MIDI  wavelength  range  from  8.0 /nn  to  1 3.0 //m  were 
used.  The  detector  masks  were  calculated  by  the  procedure  of 
MyMidiGui,  and  were  used  for  both  the  MIA  and  the  EWS  analy¬ 
sis.  The  results  obtained  by  MIA  and  EWS  were  found  to  be  vir¬ 
tually  identical  for  all  datasets.  We  chose  to  use  in  the  following 
the  results  based  on  the  EWS  analysis.  It  is  based  on  a  coherent 
integration  of  the  fringe  signal  after  numerical  compensation  of 
optical  path  differences  in  each  scan. 

Calibrated  visibility  spectra  for  S  Ori  were  derived  using  the 
transfer  function  obtained  for  our  calibrator  HD  50778.  The  scat¬ 
ter  of  all  obtained  transfer  functions  during  the  night  was  used 
to  estimate  its  uncertainty.  The  latter  depends  on  the  sample  of 
calibrations  stars  available  per  night,  and  we  assumed  a  mini¬ 
mum  uncertainty  of  the  transfer  function  of  0.1  based  on  nights 
when  many  values  were  obtained.  The  total  error  of  our  cali¬ 
brated  S  Ori  visibility  spectrum  includes  the  errors  of  the  raw 
coherence  factors  computed  by  EWS,  the  uncertainty  of  the  cali¬ 
brator  diameters,  and  the  scatter  of  the  transfer  function  during 
the  night.  The  total  error  is  dominated  by  systematic  uncertain¬ 
ties  for  the  height  of  the  interferometric  transfer  function,  so  that 
the  relative  accuracy  of  the  visibility  values  among  the  spectral 
channels  is  better  than  indicated  by  our  total  error  bars. 

We  absolutely  calibrated  the  photometric  spectrum  of  S  Ori 
for  each  data  set  using  the  interferometric  calibrator  HD  50778 

2  http : //www . ster . kuleuven . ac . be/~ti j 1/ 
MIDI_calibration/mcc . txt 

3  http://www.strw.leidenuniv.nl/~koehler/MIDI 

4  http : //www . sc . eso . org/~chummel/midi/midi . html 


also  as  the  photometric  calibration  star.  The  absolutely  cali¬ 
brated  photometric  spectrum  of  HD  50778  was  not  available 
directly.  We  used  instead  the  average  of  the  available  spectra 
of  HD  20644  and  HD  87837  (Cohen  et  al.  1999),  scaled  with 
the  IRAS  12/rm  flux  (Beichmann  et  al.  1988)  to  the  level  of 
HD  50778.  These  stars  have  the  same  spectral  type  (K4III)  and 
very  similar  effective  temperatures  as  HD  50778  (4058  K  and 
4094  K  versus  4049  K).  In  addition,  we  calibrated  the  photomet¬ 
ric  spectrum  of  HD  50778  using  the  spectro-photometric  calibra¬ 
tor  HD  31421,  which  was  observed  close  in  time  during  the  night 
2005-03-03,  and  verified  that  our  synthetic  HD  50778  spectrum 
is  valid.  We  obtained  integrated  8.0- 1 3.0  pm  S  Ori  flux  densi¬ 
ties  for  our  MIDI  epochs  A  to  D  of  187  +  12  Jy,  152  +  17  Jy, 
198  +  18  Jy,  and  203  +  13  Jy.  For  comparison,  the  integrated 
IRAS  LRS  flux  density  for  the  same  bandpass  (8.0- 1 3.0pm) 
is  169  Jy. 

The  MIDI  calibrated  visibility,  as  well  as  calibrated  flux  val¬ 
ues,  are  detailed  in  Table  A.l. 

3.2.  MIDI  data  reduction  results 

As  an  alternative  and  probably  more  intuitive  quantity  compared 
to  the  visibility  values,  we  calculated  the  equivalent  UD  diam¬ 
eters.  The  equivalent  UD  diameter  was  computed  for  each  pro¬ 
jected  baseline  length  and  spectral  channel  separately  and  de¬ 
fined  as  the  diameter  of  that  uniform  disc  that  results  in  the  same 
visibility  value  at  each  separate  data  point.  Note  that  the  true 
intensity  profile  is  not  expected  to  be  a  uniform  disc.  As  a  re¬ 
sult,  the  equivalent  UD  diameters  for  different  projected  base¬ 
line  lengths  are  not  expected  to  coincide  (not  even  for  a  cir¬ 
cularly  symmetric  intensity  profile).  However,  these  values  are 
fully  equivalent  to  the  calibrated  visibility  values  and  can  pro¬ 
vide  a  first  rough  estimate  of  the  characteristic  size  of  the  target. 

Figures  2-5  show,  for  MIDI  epochs  A-D,  respectively,  the 
obtained  MIDI  flux  spectra,  the  MIDI  visibility  values  as  a  func¬ 
tion  of  wavelength,  the  MIDI  visibility  values  as  a  function  of 
spatial  frequency  for  the  example  of  three  spectral  channels  and 
the  equivalent  UD  diameter  values.  Also  shown  are  the  synthetic 
visibility  and  flux  values  of  the  atmosphere  and  dust  models  de¬ 
scribed  below  in  Sect.  3.3  in  each  panel,  as  well  as  the  intensity 
profiles  of  these  models  for  the  example  of  three  bandpasses  in 
the  bottom  panels. 

The  general  shape  of  our  S  Ori  visibility  functions,  or  equiv¬ 
alent  UD  diameter  values,  in  Figs.  2-5  is  qualitatively  simi¬ 
lar  to  the  MIDI  data  of  the  Mira  variable  RR  Sco  described  in 
Ohnaka  et  al.  (2005).  Objects  such  as  these  are  characterized  by 
a  partially  resolved  stellar  disc,  including  atmospheric  molecu¬ 
lar  layers  that  are  optically  thick  at  mid-infrared  wavelengths, 
surrounded  by  a  spatially  resolved,  optically  thin  dust  shell. 
The  visibility  increases  between  about  8-9  /mi,  corresponding 
to  a  quasi-constant  equivalent  UD  diameter  that  is  roughly  twice 
the  continuum  photospheric  size.  Here,  the  observed  intensity 
is  dominated  by  radiation  from  optically  thick  molecular  lay¬ 
ers  at  distances  up  to  roughly  two  continuum  photospheric  radii. 
Beyond  about  9- 1 0  pm,  the  visibility  function  flattens,  corre¬ 
sponding  to  an  increasing  equivalent  UD  diameter.  In  this  wave¬ 
length  region,  spatially  resolved,  optically  thin  radiation  from 
the  dust  shell  starts  to  represent  a  considerable  part  of  the  total 
intensity.  In  addition,  extinction  of  the  stellar  light  by  the  dust 
shell  becomes  important.  At  the  longest  wavelengths  of  MIDI, 
spatially  resolved  radiation  from  the  dust  shell  dominates  the 
measured  intensity.  A  detailed  model  of  S  Ori’s  atmosphere  and 
dust  shell  is  described  in  Sects.  3.3  and  3.4. 
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Fig.  2.  Results  of  our  VLTI/MIDI  8- 13 /mi  interferometry  of 
S  Ori  at  epoch  A,  stellar  phase  0.42.  Shown  are  our  MIDI  data 
compared  to  the  best-fitting  model  with  parameters  listed  in 
Table  3.  In  each  panel,  the  x-symbols  with  error  bars  denote 
the  measured  values,  the  solid  lines  the  synthetic  values  of  the 
overall  model  consisting  of  the  stellar  model  atmosphere  and 
the  dust  shell.  The  dotted  lines  denote  the  contribution  from 
the  attenuated  atmosphere  model  alone,  and  the  dashed  lines 
the  contribution  of  only  the  dust  shell  (attenuated  input  radia¬ 
tion  subtracted),  a)  The  calibrated  MIDI  flux  spectrum,  b)  The 
calibrated  MIDI  visibility  values  as  a  function  of  wavelength, 
c)  The  calibrated  MIDI  visibility  values  as  a  function  of  spatial 
frequency  for  the  example  of  three  spectral  channels  at  wave¬ 
lengths  8.5 /rm,  10.5 /i m,  and  12.5  // m.  d)  The  equivalent  UD  di¬ 
ameter  values.  Also  indicated  are  the  photospheric  angular  di¬ 
ameter  ©phot  and  the  inner  dust  shell  boundary  0in  as  obtained 
by  our  model  fit,  as  well  as  the  mean  42.8  GHz  and  43.1  GHz 
SiO  ring  angular  diameters  (Sect.  4).  e)  The  model  intensity 
profiles  at  spectral  bands  8-9 /rm,  10-11/rm,  and  12-13/rm. 
The  bottom  abscissa  shows  the  radius  as  a  function  of  the  pho¬ 
tospheric  radius  Rphot,  and  the  top  abscissa  the  corresponding 
angular  scale  obtained  from  our  model  fit.  The  arrows  indicate 
the  location  of  the  inner  boundary  radius  of  the  dust  shell  Rm, 
as  well  as  of  the  mean  43.1  GHz  and  42.8  GHz  SiO  maser  spot 
radii  R43.1  and  ^42.8- 


3.3.  MIDI  data  modeling 

Mid-infrared  photometric  and  interferometric  data  of  Mira  stars 
are  in  general  sensitive  to  details  of  the  structure  of  the  stellar 


atmosphere  in  which  molecular  layers  of  various  geometrical 
and  optical  characteristics  (e.g.,  Tej  et  al.  2003a)  lie  above 
the  continuum-forming  photosphere,  as  well  as  to  the  sur¬ 
rounding  dust  shell  (cf.  e.g.,  Ohnaka  et  al.  2005,  2006b; 
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Fig.  3.  Results  of  our  VLTI/MIDI  8-13/rm  interferometry  of  S  Ori  at  Fig.  4.  Results  of  our  VLTI/MIDI  8-13  /rm  interferometry  of  S  Ori  at 
epoch  B,  stellar  phase  0.55.  As  Fig.  2.  epoch  C,  stellar  phase  1.16.  As  Fig.  2. 
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Fig.  5.  Results  of  our  VLTI/MIDI  8-13/rm  interferometry  of  S  Ori  at 
epoch  D,  stellar  phase  1.27.  As  Fig.  2. 


Ireland  &  Scholz  2006;  Weiner  et  al.  2006).  We  used  the  P  and 
M  model  atmosphere  series  (Ireland  et  al.  2004a,b)  to  describe 


the  stellar  atmosphere  including  the  continuum  photosphere 
and  overlying  molecular  layers,  and  added  an  ad-hoc  radiative 
transfer  model  of  the  dust  shell.  The  details  of  our  modeling  are 
described  in  the  following. 

3.3.1 .  Atmosphere  model  for  S  Ori 

Few  dynamic  atmosphere  models  for  oxygen-rich  Mira  stars 
are  available  that  include  the  effects  of  molecular  layers  (see, 
e.g.,  Scholz  2001,  2003).  The  P  and  M  model  series  (Ireland 
et  al.  2004a, b)  are  complete  self-excited  dynamic  model  atmo¬ 
sphere  of  Mira  stars,  and  they  have  successfully  explained  a  lot 
of  observational  data  (e.g.;  Hofmann  et  al.  2000,  2001,  2002; 
Woodruff  et  al.  2004;  Fedele  et  al.  2005).  These  model  series 
have  been  constructed  to  match  the  prototype  oxygen-rich  Mira 
stars  o  Cet  and  R  Leo.  The  P  and  M  series  differ  with  respect  to 
the  mass  of  the  so-called  “parent  star”,  which  is  the  hypotheti¬ 
cal  non-pulsating  equivalent  of  the  pulsating  Mira  variable.  The 
geometric  pulsation  of  the  Mira  occurs  around  the  parent  star’s 
Rosseland  radius  R  =  R(tross  =  1)  (Ireland  et  al.  2004a,b).  The 
parent  star  has  solar  metallicity,  luminosity  L/L®  =  3470,  pe¬ 
riod  332  days,  mass  M/MQ  =  1.0  (P  series)  and  1.2  (M),  radius 
Rp/Rq  =  241  (P)  and  260  (M)  (Hofmann  et  al.  1998).  The  mod¬ 
erately  larger  radius  of  the  M-series  parent  star  leads  to  a  slightly 
lower  effective  temperature  for  the  parent  star  and  systematically 
lower  phase-dependent  effective  temperatures  for  the  pulsating 
Mira  compared  to  the  P  series.  The  M  models  tend  to  exhibit 
less  pronounced  cycle-to-cycle  variations  than  the  P  models  and 
to  have  more  compact  atmospheres  (Ireland  et  al.  2004b).  Effects 
on  interferometric  diameter  measurements  resulting  from  the 
differences  between  the  P  and  M  model  series  are  subtle.  When 
comparing  P  and  M  models  of  similar  phases  to  near-infrared 
interferometric  data  of  R  Leo,  the  resulting  Rosseland  angular 
diameter  values  were  found  to  agree  within  $1.5%  (Fedele  et  al. 
2005). 

Compared  to  o  Cet  and  R  Leo,  for  which  the  P  and  M  model 
series  are  designed,  S  Ori  is  a  slightly  cooler  Mira  variable 
(M6.5-M9.5  versus  M5-M9/M6-M9.5,  Samus  et  al.  2004) 
with  a  longer  period  (420d  versus  332d/310d,  Samus  et  al. 

2004) ,  a  larger  main  sequence  precursor  mass  (1.3  MQ  ver¬ 
sus  1.0Mo/1.2Mo,  Wyatt  &  Cahn  1983),  and  a  larger  radius 
(mean  continuum  photospheric  radii  roughly  R  ~  450  RQ  ver¬ 
sus  ~350Ro;  based  on  BW05,  Woodruff  et  al.  2004;  Fedele  et  al. 

2005) .  However,  when  scaled  to  variability  phases  between  0 
and  1  and  to  the  corresponding  angular  size  on  sky,  the  general 
model  results  are  not  expected  to  be  dramatically  different  for 
S  Ori  compared  to  Mira  stars,  such  as  o  Cet  and  R  Leo.  Hereby, 
the  M  series  is  expected  to  be  better-suited  to  describing  S  Ori 
than  the  P  series  because  of  its  higher  parent-star  mass.  We  con¬ 
cluded  that  the  use  of  the  M  series  is  currently  the  best  available 
option  for  describing  S  Ori’s  atmosphere  including  the  contin¬ 
uum  photosphere  and  overlying  atmospheric  molecular  layers. 

3.3.2.  Dust-shell  model  forS  Ori 

Except  for  the  dust-enshrouded  very-long-period  Mira  IW  Hya 
(Jeong  et  al.  2003),  models  that  self-consistently  combine 
detailed  dynamic  atmosphere  models  and  the  formation  of 
dust  shells  are  currently  not  readily  available  for  oxygen-rich 
AGB  stars.  However,  the  first  attempts  in  this  direction  have  just 
been  made  (Ireland  &  Scholz  2006,  Woitke  2006). 

Dust  shells  around  AGB  stars  have  often  been  modeled 
using  ad-hoc  radiative  transfer  calculations  (see;  for  instance; 
Dane  hi  et  al.  1994;  Lorenz -Martins  &  Pompeia  2000;  Ohnaka 
et  al.  2005,  2006a),  and  we  followed  this  approach  as  the  only 
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currently  available  option  for  describing  S  Ori’s  dust  shell.  We 
employed  the  Monte  Carlo  radiative  transfer  code  mcsim_mpi 
(Ohnaka  et  al.  2006a)  to  calculate  the  dust  shell  model  of  S  Ori. 
For  the  central  radiation  source  of  the  radiative  transfer  model, 
we  used  the  spectral  energy  distribution  (SED)  of  the  respective 
M  model  atmosphere  for  the  wavelength  range  from  0.32 //m 
to  23 //m  and  extended  it  by  a  blackbody  approximation  for 
longer  wavelengths.  Dust  temperatures  were  calculated  in  radia¬ 
tive  equilibrium. 

Our  MIDI  spectra  of  S  Ori  do  not  exhibit  any  prominent 
silicate  feature,  indicating  that  silicate  grains  are  not  a  major 
constituent  of  the  dust  shell.  Lorenz-Martins  &  Pompeia  (2000) 
have  investigated  the  envelopes  of  3 1  oxygen-rich  AGB  stars  us¬ 
ing  the  IRAS  LRS  spectra.  They  classify  the  objects  into  three 
groups,  whose  IRAS  LRS  spectra  can  be  modeled  using  either 
silicate  grains,  AFO?  grains,  or  a  mix  thereof.  S  Ori  is  part  of 
their  study  and  belongs  to  the  group  of  objects  that  can  be  mod¬ 
eled  well  with  AI2O3  grains  alone.  They  show  that  a  blend  of 
alumina  AFO3  grains  (Koike  et  al.  1995)  for  A  <  8  pm  and 
porous  amorphous  AFO3  grains  (Begemann  et  al.  1997)  for 
A  >  8  //m  gives  the  best  agreement  with  the  IRAS  spectra,  so 
we  followed  this  approach  (as  in  Ohnaka  et  al.  2005  for  the  case 
of  RR  Sco).  Likewise,  we  considered  small  amounts  of  warm 
silicates  (Ossenkopf  et  al.  1992)  following  the  same  works.  The 
grain  size  was  set  to  0.1  p m  for  both  AFO3  and  silicates.  We 
used  the  same  inner  boundary  radius  and  density  distribution  for 
both  dust  species.  The  amount  of  dust  of  each  species  was  char¬ 
acterized  by  the  optical  depths  at  a  wavelength  of  0.55  //m,  de¬ 
noted  rpIAFOr)  and  Ty(silicate).  The  density  distribution  p(r) 
was  described  by  a  single  power  law  with  index  p  (p(r )  oc  r~p). 
We  used  the  inner  boundary  radii  in  units  of  the  stellar  photo- 
spheric  radius  All  radiative  transfer  models  of  the  dust 

shell  were  calculated  with  an  outer  radius,  so  that  the  shell  thick¬ 
ness  is  R0ut/Rin  =  1000. 


3.3.3.  The  main  model  parameters 

The  stellar  atmospheric  structure  including  the  effects  from  at¬ 
mospheric  molecular  layers  is  fully  described  by  the  choice 
of  the  specific  model  of  the  M  series.  The  M  series  consists 
of  20  models  for  different  phase  and  cycle  combinations  as  de¬ 
tailed  in  Ireland  et  al.  (2004b).  Here,  we  used  a  grid  of  the 
six  M  models  of  cycle  2  of  the  series,  which  are  M20  (model 
visual  phase  0.05),  M21n  (0.10),  M22  (0.25),  M23n  (0.30), 
M24n  (0.40),  and  M25  (0.50).  These  models  cover  the  range  of 
our  observations  at  phases  0.16  (epoch  C),  0.27  (epoch  D),  0.42 
(epoch  A),  and  0.55  (epoch  D).  The  absolute  model  phases  with 
respect  to  observations  are  uncertain  by  at  least  0.1,  with  rela¬ 
tive  uncertainties  of  0.01  to  0.02  between  different  models  of  the 
series  (Ireland  et  al.  2004b).  Likewise,  the  absolute  assignment 
of  visual  phases  to  our  epochs  of  observation  is  uncertain  by 
about  0.1  (see  Sect.  3).  Lor  each  of  these  6  atmosphere  models, 
whose  SEDs  are  used  as  central  source,  we  calculated  a  grid  of 
radiative  transfer  models  of  the  dust  shell.  These  grids  included 
all  combinations  of  optical  depths  rylAFCF)  =  1.0, 1.5,  2.0, 2.5, 
3.0;  Ty(silicate)  =  0%,  10%  of  rylAFCF);  inner  boundary  radii 
Rin/Rphot  =  1.8,  2.0,  2.2,  2.4,  2.6;  and  gradients  p  =  2.0,  2.5, 
3.0,  3.5,  4.0. 

Any  overall  model  was  thus  defined  with  this  set  of  five 
independent  parameters:  (1)  choice  of  the  dust-free  M  model, 
(2)  TyfAFCF),  (3)  ry(silicate),  (4)  Rm,  (5)  p.  We  used  the  well- 
defined  continuum  photospheric  radius  at  A  =  1.04//m  R1.04 
of  the  respective  dust-free  M  model  to  characterize  the  photo- 
spheric  radius  Rphod  he.  we  set  Rphot  =  Ri.04-  Any  overall  model 
is  rigid;  i.e.  the  distance  of  any  model  layer  from  Rphot  is  fixed 


with  the  parameters  given  above.  The  corresponding  angular  di¬ 
ameter  ©phot  was  treated  as  the  only  fit  parameter  when  compar¬ 
ing  any  such  overall  model  to  our  MIDI  data. 


3.3.4.  Synthetic  visibility  and  flux  values 
Stellar  contribution  Monochromatic  CLVs  of  the  M  models 
were  calculated  at  3000  equally-spaced  frequencies  for  the 
wavelength  range  8  pm  <  A  <  1 3  pm  as  described  in  Ohnaka 
et  al.  (2006b),  covering  radii  r  between  0  and  Rj'j'j  =  5  x  Rp. 
These  CLVs  were  calculated  with  H2O  and  SiO  (28SiO,  29SiO, 
and  30SiO)  lines  included.  The  line  lists  used  for  the  calculations 
were  the  same  as  those  used  for  modeling  RR  Sco  by  Ohnaka 
et  al.  (2005);  and  as  in  the  case  of  RR  Sco,  we  adopted  silicon 
isotope  ratios  of  28Si/29Si  =  15  and  28Si/30Si  =  20.  The  H20 
and  SiO  line  opacities  were  computed  assuming  local  thermo¬ 
dynamical  equilibrium  (LTE)  and  a  Voigt  profile  with  a  micro- 
turbulent  velocity  of  3  km  s_1 . 

The  visibility  values  based  on  the  stellar  contribution  were 
calculated  for  each  spectral  channel  c  of  MIDI  as 


ystar 
r  c 


rVmaxF)  rSom  rstar/  \  p-t(v) 

Jvmin(c)  Jo  A 


JQ(n  (B/A))  r  dr  dv 


rVmax(c) 
-^min  (c) 


e~T(v>r  dr  dv 


(1) 


Here,  B  is  the  projected  baseline  length  used  for  the  observa¬ 
tion.  Note  that  this  computation  is  valid  for  the  case  of  the 
EWS  data  reduction,  which  coherently  integrates  the  fringe  sig¬ 
nal  to  directly  obtain  the  visibility  modulus.  Lor  the  MIA  analy¬ 
sis,  which  is  a  power  spectrum  analysis  integrating  the  squared 
visibility  modulus,  the  integration  over  frequency  would  need  to 
be  performed  over  the  squared  monochromatic  visibility  values 
(cf.  Wittkowski  et  al.  2006).  The  observed  stellar  flux  contribu¬ 
tion  /cstar  was  computed  for  each  of  MIDI’s  spectral  channels  as 


/•star  _  o  _ 
Jc  ~2rr 


/^\star 

wout 


n/2 


206265 


/''Vmax(c)  p 
Jymin(c)  JO 


Iy  M(r)  e“T(vV  dr  dv.  (2) 


The  spectral  channels  of  MIDI  were  assumed  to  have  a  rectangu¬ 
lar  sensitivity  function  between  frequencies  vm;n  and  vmax.  The 
term  e  'T<v>  describes  the  extinction  by  the  dust  shell.  The  opti¬ 
cal  depth  t(v)  was  calculated  based  on  tabulated  absorption  and 
scattering  coefficients  (2abs  and  Gsca  of  the  two  dust  species  i  as 


t(v)  =  ^  r'(vo) 

i 


<2Lbs(y)  +  QL(y) 

6a  bs(vo)  +  GscaOo)' 


(3) 


where  we  used  vq  corresponding  to  Ag  =  0.55  pm. 


Dust  contribution  Using  the  dust  temperatures  and  monochro¬ 
matic  mean  intensities  computed  by  a  Monte  Carlo  run  with 
mcsim_mpi  for  a  given  set  of  dust  shell  parameters,  we  cal¬ 
culated  monochromatic  CLVs  of  the  dust  shell  (excluding  the 
contribution  of  the  central  star)  with  the  ray-tracing  method  as 
described  in  Ohnaka  et  al.  (2006a,  2007).  The  radii  r  covered 
the  range  between  0  and  25  Rphot,  and  we  denote  the  outer¬ 
most  model  radius,  where  we  set  the  intensity  to  0,  by  R^;st  = 
25  x  Rphot-  The  outer  radius  of  25  Rphot  was  chosen  because  the 
MIDI  interferometric  field  of  view  of  ~250  mas  corresponds  to 
25  times  the  expected  photospheric  angular  diameter  of  S  Ori 
of  ~  10  mas.  We  verified  that  changing  the  outer  model  radius 
from  25  Rphot  to  50  Rphot  changed  the  diameter  results  by  less 
than  0.3%.  This  indicates  that  the  contribution  from  dust  to  the 
overall  CLV  outside  of  25  photospheric  radii  can  be  neglected 
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Table  3.  Fit  results  of  our  MIDI  measurements  of  S  Ori. 


Epoch  (®vis) 

A  (0.42) 

B  (0.55) 

C  (1.16) 

D  (1.27) 

Best  model 

M22 

M24n 

M23n 

M21n 

®Model 

0.25 

0.40 

0.30 

0.10 

0Model  —  0Vis 

-0.17 

-0.15 

+0.14 

-0.17 

Best  Rin/Rphot 

1.8 

2.0 

2.2 

2.4 

Best  /HALCb) 

2.5 

2.5 

1.5 

1.5 

Best  Tp(Silicate) 

0 

0 

0 

0 

Best  p 

3.5 

3.5 

3.0 

2.5 

©Phot  [mas] 

9.00 

7.92 

9.66 

9.48 

xl 

0.16 

0.10 

0.43 

0.35 

n 

62 

186 

124 

62 

^"formal  [mas] 

0.10 

0.05 

0.02 

0.07 

^"model  [mas] 

0.30 

0.09 

0.06 

0.34 

and  that  the  exact  choice  of  the  outermost  model  radius  is  not 
critical  for  our  calculation.  Visibility  (Vdust)  and  flux  (/cdust)  val¬ 
ues  for  each  spectral  channel  of  MIDI  were  computed  in  the 
same  way  as  for  the  dust-free  model  atmosphere,  but  without 
the  extinction  term  e_T(v). 

Comparison  of  visibility  and  flux  values  to  our  data  Finally, 
the  flux  and  visibility  values  of  the  overall  model  were 
computed  as 

/total=/,tar+/dust  (4) 

and 

rstar  rdust 

T/total  •' c  x/star  ,  J c  x/dust  /r\ 

V c  ~  /-total  V c  +  /-total  Yc  '  V-V 

Jc  J  c 

This  simple  addition  of  the  CLVs  of  the  dust-free  model  atmo¬ 
sphere  of  the  M  series  and  of  the  pure  dust-shell  (i.e.  after  sub¬ 
traction  of  the  contribution  by  the  attenuated  input  radiation)  is 
fully  valid  as  long  as  the  stellar  atmosphere  and  dust  shell  are 
spatially  separated.  In  fact,  there  may  be  some  overlap  of  the  ex¬ 
tended  low-intensity  atmospheric  molecular  layers  of  the  dust- 
free  atmosphere  model  and  the  dust  shell  in  our  overall  mod¬ 
els.  However,  we  verified  that  the  inner  boundary  radius  of  the 
dust  shell  is  clearly  larger  than  that  of  the  layer  at  which  the  at¬ 
mosphere  becomes  optically  thick  at  mid-infrared  wavelengths. 
Ireland  &  Scholz  (2006)  find  for  similar  conditions  that  the  gas 
temperatures  for  a  dust-free  M  model  and  an  M  model  that  in¬ 
cludes  re-heating  of  gas  by  dust,  differ  by  less  than  50  K  for  radii 
up  to  about  3  parent-star  radii. 

For  any  given  overall  model  separately,  we  optimized  the  an¬ 
gular  diameter  ©phot  using  a  standard  least-square  fit  that  mini¬ 
mizes  the^f2  value  between  synthetic  and  observed  visibility  and 
flux  values.  Hereby,  ©phot  was  treated  as  the  only  free  parameter. 
The  resulting  reduced^2  values  among  our  full  grid  of  different 
overall  models  were  compared  in  order  to  find  the  best-fitting  set 
of  model  parameters.  This  analysis  was  separately  performed  for 
each  of  our  epochs  A-D. 


3.4.  MIDI  model  results 

Table  3  lists  for  each  of  our  MIDI  epochs  the  parameter  set 
of  the  best-fitting  overall  model  of  our  grid.  Listed  are  for 
each  of  our  4  epochs  the  best-fitting  atmosphere  model  of  the 
M  series;  its  model  phase;  the  difference  between  observed 
and  model  phase;  the  best-fitting  parameters  of  the  dust  shell 
(tvIALOs),  Ty(Silicate),  ry(Silicate),  p)\  the  best-fitting  angular 
diameter  0phofi  the  reduced^2  value;  the  number  of  data  points 
used  (n);  and  the  formal  and  model  errors.  The  model  errors  are 


Table  4.  Reduced^2  values  obtained  for  our  grid  of  overall  models. 


Epoch 
M  model 


Rin/Rphot 


Ty(Al203) 


Ty(Silicate)/ 

TyTALCF) 

P 


M20 


A 

0.27 


B 

0.16 


C 

0.52 


D 

0.38 


M2  In 
M22 
M23n 
M24n 
M25n 
1.8 
2.0 
2.2 
2.4 
2.6 
1.0 


0.26 

0.16 

0.17 

0.20 

0.19 

0.16 

0.19 

0.21 

0.24 

0.29 

/ 


0.14 

0.17 

0.12 

0.10 

0,17 

0.10 

0.10 

0.11 

0.12 

0.18 

/ 


0.53 

0.81 

0.43 

0.66 

0.53 

0.56 

0.52 

0.43 

0.50 

0.57 

0.77 


0.35  | 
0.51 
0.61 
1.04 
0.58 
0.54 
0.46 
0.39 
0.35 
0.38 
0.48 


1.5 
2.0 

2.5 
3.0 
0% 

10% 

2.0 


0.18 

0.17 

0.16 

0.18 

0.16 

0.19 

0.31 


0.11 

0.11 

0.10 

0.11 

0.10 

0.12 

0.20 


0.43  | 
0.56 
0.80 

/ 

0.43 

0.52 

0.87 


0.35  | 
0.54 
0.71 

/ 

0.35 

0.38 

0.38 


2.5 
3.0 

3.5 
4.0 


0.20 

0.17 

0.16 

0.20 


0.14 

0.10 

0.10 

0.12 


0.54 

0.43 

0.50 

/ 


0.35  | 
0.39 
0.51 

/ 


derived  as  the  standard  deviations  of  ©phot  values  based  on  differ¬ 
ent  overall  models  of  our  grid  that  are  consistent  at  the  1  cr  level 
based  on  application  of  the  F-test.  The  model  errors  of  epochs  A 
and  D  are  larger  than  those  of  epochs  B  and  C,  which  can  be  ex¬ 
plained  by  the  use  of  only  one  baseline  configuration  at  these 
epochs.  Table  4  lists  for  each  epoch  and  for  each  grid  parameter 
the  best  xl  value  that  could  be  obtained  by  any  combination  of 
all  other  parameters.  The  best  xl  value  for  each  epoch  is  marked 
by  boxes. 

The  synthetic  visibility  and  flux  values  corresponding  to  the 
best-fitting  overall  models  are  shown  in  Figs.  2-5  in  comparison 
to  the  measured  values.  Panels  (d),  which  show  the  equivalent 
UD  diameter  values,  also  indicate  the  photospheric  angular  di¬ 
ameter  0phot  and  the  inner  dust  shell  boundary  0jn  derived  from 
the  model  fit.  Panels  (e)  show  the  intensity  profiles  of  the  over¬ 
all  models  that  fit  our  MIDI  data  best  for  three  bandpasses  at 
8-9/rm,  10-lljum,  and  12-13 pm.  These  panels  have  a  bottom 
and  a  top  abscissa  indicating  (bottom)  the  radius  as  a  function 
of  the  model  photospheric  radius  Rphot  and  (top)  the  correspond¬ 
ing  angular  scale  using  the  best-fitting  angular  diameter  ©phot- 
Panels  (e)  also  indicate  the  inner  dust  shell  radius  Rm,  as  well 
as  the  mean  SiO  maser  ring  radii  R42.8  and  R43.1,  which  are  de¬ 
rived  in  Sect.  4.  A  two-dimensional  pseudo-color  image  of  the 
mid-infrared  model  intensity  showing  the  molecular  layers  and 
the  dust  shell  is  displayed  in  Fig.  7.  The  photospheric  disk  is 
represented  by  a  light  blue  color.  Overlaid  are  the  images  of  the 
SiO  maser  radiation  described  in  Sect.  4. 

The  photospheric  diameter  ©Phot  The  obtained  photospheric 
diameter  values  ©phot  at  phases  0.42,  0.55,  1.16,  and  1.27  are 
9.0mas  +  0.3 mas,  7.9 mas  +  0.1  mas,  9.6mas  +  0.1  mas,  and 
9.5  mas  ±  0.4  mas,  respectively.  Here,  the  error  bars  include  the 
formal,  as  well  as  the  model,  errors  from  Table  3.  Our  measure¬ 
ments  indicate  significant  changes  of  the  photospheric  diame¬ 
ter  Ophot  as  a  function  of  stellar  phase.  A  detailed  comparison 
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of  ©phot  to  the  inner  dust  shell  radius  and  the  SiO  maser  ring 
radii  follows  in  Sect.  5. 

In  BW05,  we  estimated  S  Ori’s  photospheric  angular  diam¬ 
eter  at  phase  0.73  to  ~9.2  mas,  based  on  our  /if -band  UD  diame¬ 
ter  measurements  and  using  correction  factors  from  UD  to  pho¬ 
tospheric  diameters  from  Ireland  et  al.  (2004a,b).  This  value  is 
consistent  with  the  photospheric  diameters  derived  here,  as  the 
estimate  of  9.2  mas  at  phase  0.73  from  BW05  lies  as  expected 
in  between  the  values  of  7.9  mas  at  phase  0.54  and  9.6  mas  at 
phase  1.16  derived  in  the  present  paper. 

The  K -band  UD  diameter  of  10.54 mas  +  0.68  mas  by 
van  Belle  et  al.  (1996)  was  obtained  at  phase  0.56.  At  phases 
close  to  this  value,  the  M  series  results  in  a  ratio  between  the 
K -band  UD  diameter  and  the  photospheric  diameter  of  —  1 .35 
(Fig.  7  of  Ireland  et  al.  2004b).  With  this  correction,  the 
photospheric  diameter  corresponding  to  this  measurement  at 
phase  0.56  would  be  7.8  mas  ±  0.5  mas,  and  would  thus  be  fully 
consistent  with  our  value  of  7.9  mas  +  0.1  mas  at  the  same  phase 
(Epoch  B  at  phase  0.54). 

Millan-Gabet  et  al.  (2005)  obtained  a  /('-hand  UD  diameter 
of  9.6  mas  +  0.2  mas  at  phase  0.1.  At  this  phase,  the  ratio  be¬ 
tween  //-band  UD  diameter  and  photospheric  diameter  is  pre¬ 
dicted  to  be  close  to  unity  (Ireland  et  al.  2004b).  This  measure¬ 
ment  thus  also  agrees  well  with  our  photospheric  diameter  of 
9.6  mas  ±0.4  mas  at  phase  1.16  (epoch  C). 

The  good  agreement  of  the  photospheric  angular  diameters 
derived  from  our  MIDI  measurements  in  the  present  paper  with 
previous  near-infrared  diameter  measurements  increases  confi¬ 
dence  in  the  validity  of  our  modeling  using  the  M  model  series 
complemented  by  a  radiative  transfer  model  of  the  dust  shell. 

The  choice  of  the  M  model  The  intensity  profile  of  the  ex¬ 
tended  atmosphere,  and  thereby  the  choice  of  the  best-fitting 
model  out  of  the  M  series,  is  mainly  constrained  by  the  shape  of 
the  visibility  functions,  in  particular  for  longer  baseline  lengths, 
in  the  wavelength  range  ~8— 9  fin t,  where  the  contribution  by  the 
dust  shell  is  relatively  low. 

The  differences  in  visual  stellar  phase  of  the  best-fitting 
M  model  compared  to  the  visual  phase  at  the  epochs  of  obser¬ 
vation  are  -0.17,  -0.15,  +0.14,  and  -0.17,  respectively.  These 
differences  may  be  due  to  the  uncertainties  of  the  absolute  as¬ 
signments  of  visual  phases  to  the  model  series  as  well  as  to  the 
dates  of  observation,  which  are  both  on  the  order  of  0.1.  The 
relative  phase  values  are  expected  to  agree  better,  which  is  true 
for  epochs  A,  B,  and  D.  Differences  between  these  phases  might 
also  occur  owing  to  the  different  radius  and  pulsation  period  of 
the  model  series’  parent  star  compared  to  S  Ori.  These  parame¬ 
ters  might  have  an  effect  on  the  extension  of  the  atmosphere,  so 
that  the  best-fitting  extension  of  the  atmosphere  might  belong  to 
a  model  at  a  different  visual  phase.  An  effect  such  as  this  might 
explain  the  clearly  different  phase  difference  between  observa¬ 
tion  and  model  at  epoch  C  compared  to  the  three  other  epochs. 

The  good  agreement  of  the  shape  of  the  synthetic  visi¬ 
bility  function  with  our  observations  at  the  wavelength  range 
from  8  fim  to  9  fim  and  the  simultaneous  agreement  with  the 
measured  flux  level  again  confirms  that  the  M  model  series  is 
usable  for  S  Ori. 

As  an  indication  of  the  extension  of  the  molecular  layers 
of  the  best-fitting  M  models,  the  50%/10%  intensity  radii  of 
the  A-band  (8-13 pm)  model  CLVs  are  1.5/1.6,  1.8/1.9,  1.6/1, 7, 
and  1.4/1.9  photospheric  radii //phot  for  epochs  A-D,  i.e.  for  mod¬ 
els  M22,  M24n,  M23n,  and  M2  In,  respectively. 


The  dust  shell  chemistry  The  ratio  of  the  number  of 
AFO(  grains  to  that  of  silicate  grains  is  mostly  constrained  by 
the  shape  of  the  flux  and  visibility  spectra  at  the  silicate  fea¬ 
ture  close  to  1 0  fim.  We  obtain  best  fits  to  our  S  Ori  data  with  a 
dust  shell  only  consisting  of  AFO  i  grains.  The  addition  of  small 
amounts  of  silicate  grains  does  not  improve  the  fits  for  any  of  our 
epochs.  While  our  grid  in  Table  4  only  includes  silicate  grains 
with  the  same  inner  boundary  radius  as  for  the  AFO?  grains,  we 
also  compared  to  our  data  the  best-fitting  models  from  Table  3 
with  silicate  grains  added  at  larger  inner  boundary  radii  of  4 
and  6  photospheric  radii.  These  models  lead  to  clearly  increased 
xl  values  as  well.  In  summary,  our  data  indicate  that  silicon  is 
not  bound  in  dust,  neither  at  very  close  distances  of  ~2  Rph0t  nor 
at  larger  distances  of  4-6  //phoi-  A  dust  shell  only  consisting  of 
AFOi  grains  is  consistent  with  the  modeling  of  the  IRAS  LRS 
spectra  of  S  Ori  by  Lorenz-Martins  &  Pompeia  (2000). 

The  radial  density  profile  of  the  dust  shell  The  total  optical 
depth  of  the  dust  shell  is  mainly  constrained  by  the  height  and 
shape  of  the  MIDI  total  flux  spectrum.  The  radial  density  profile 
described  by  the  inner  boundary  radius  and  the  density  gradient 
is  mostly  constrained  by  the  shape  of  the  visibility  function  ver¬ 
sus  wavelength.  Here,  the  values  at  shorter  wavelengths,  prob¬ 
ing  warmer  dust,  constrain  the  inner  boundary  radius  best.  The 
density  gradient  is  mostly  constrained  by  the  visibility  values  at 
longer  wavelengths  and  shorter  baseline  lengths,  with  which  the 
dust  shell  is  less  resolved. 

We  obtain  best-fitting  inner  boundary  radii  of  1.8,  2.0,  2.2, 
and  2.4  photospheric  radii  Rphot  at  our  phases  0.42,  0.55,  1.16, 
and  1 .27,  respectively.  These  radii  are  well-constrained  and  do 
not  significantly  change  for  variations  in  other  model  parame¬ 
ters.  We  estimate  the  uncertainty  of  this  parameter  to  about  0.2  x 
Rphot-  The  formation  of  dust  consisting  of  AFOi  grains  at  these 
short  distances  from  the  stellar  surface  is  consistent  with  the  em¬ 
pirical  results  by  Lorenz-Martins  &  Pompeia  (2000),  as  well 
as  with  the  recent  theoretical  calculations  by  Ireland  &  Scholz 
(2006)  and  Woitke  (2006). 

Best-fitting  density  gradients  p  at  the  same  phases  are  3.5, 
3.5,  3.0,  and  2.5,  respectively,  and  best-fitting  optical  depths  ry 
are  2.5,  2.5,  1.5,  and  1.5.  Taking  the  inner  boundary  radii,  the 
density  gradients,  and  the  optical  depth  values  together,  our  fits 
indicate  more  compact  dust  shells  with  larger  optical  depth  near 
the  stellar  minimum  and  more  extended  dust  shells  with  lower 
optical  depth  after  stellar  maximum. 

Wind  models  We  used  the  radiation  transport  code  DUSTY 
(Ivezic  &  Elitzur  1997;  Ivezic  et  al.  1999)  with  its  option  of 
modeling  radiatively  driven  winds  to  compute  the  density  struc¬ 
ture  of  S  Ori’s  envelope  that  would  result  from  a  wind  model. 
We  computed  wind  models  for  each  of  our  epochs  A  to  D  us¬ 
ing  a  blackbody  with  effective  temperature  of  the  best-fitting 
M  model  as  the  central  radiation  source  and  with  the  Rm  and 
ry  values  from  Table  3.  We  used  the  same  optical  properties 
of  AFO;?  as  above,  and  we  set  Rout/Rm  to  103.  Figure  6  shows 
the  resulting  normalized  density  profiles  of  the  wind  models  to¬ 
gether  with  the  p  oc  r  1  and  p  oc  r  3  curves.  The  wind  models 
result  in  density  structures  p  oc  r-2  for  r  <;  lORphot-  Closer  to 
the  photosphere,  the  gradients  of  the  wind  density  profiles  be¬ 
come  steeper  than  p  oc  r~3.  As  our  MIDI  visibility  measure¬ 
ments  are  mostly  sensitive  to  the  inner  dust  region,  the  relatively 
high  density  gradients  given  in  Table  3  could  be  explained  by 
such  a  wind  model.  The  differences  among  the  empirically  found 
density  gradients  for  the  different  epochs  could  be  explained  by 
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Fig.  6.  Density  profiles  of  radiatively  driven  wind  models  calculated 
with  the  DUSTY  code  (Ivezic  &  Elitzur  1997;  Ivezic  et  al.  1999)  for  input 
parameters  from  Table  3,  together  with p  oc  r~2  and  p  oc  r~3  curves. 


different  mean  distances  from  which  most  mid-infrared  photons 
escape.  The  wind  models  shown  correspond  to  mass-loss  rates 
between  ~2.4  x  1 0  6  M0/yr  (Epoch  C)  and  ~3.5  x  10~6  M0/yr 
(Epoch  B).  The  indication  of  more  compact  dust  shells  with 
larger  optical  depth  near  stellar  minimum  and  more  extended 
dust  shells  with  lower  optical  depth  after  stellar  maximum  thus 
suggests  a  wind  with  increased  mass  loss  and  dust  formation 
occurring  relatively  close  to  the  stellar  surface  near  minimum 
visual  phase  and  with  an  expanded  dust  shell  after  maximum 
phase. 

Symmetry  of  the  molecular  and  dust  shells  Our  MIDI  data  at 
epochs  A  and  D  are  each  obtained  at  one  single  position  angle 
on  sky.  Epoch  B  spans  a  range  in  position  angles  of  ~40°  and 
epoch  C  of  ~20°.  Within  these  probed  ranges,  our  S  Ori  data 
can  be  described  well  by  spherically  symmetric  distributions  of 
molecular  layers  and  dust  shells  without  the  need  to  introduce 
deviations  from  circular  symmetry  to  our  modeling.  However, 
we  cannot  rule  out  asymmetries  at  wider  ranges  of  position  an¬ 
gles,  in  which  case  the  resulting  parameters  in  Table  3  are  only 
valid  along  the  position  angles  on  the  sky  given  in  Table  1 .  Note 
that  the  distributions  of  the  SiO  maser  spots  toward  S  Ori  shown 
in  BW05,  as  well  as  below  in  Sect.  4,  do  not  indicate  a  globally 
asymmetric  gas  distribution. 

4.  VLBA  measurements 

Contemporaneous  with  the  VLTI/MIDI  observa¬ 
tions,  we  observed  the  u  =  1,/  =  1-0  (43.1GHz) 
and  v  -  2,  /  =  1-0  (42.8  GHz)  SiO  maser  emission  associ¬ 
ated  with  S  Ori  (a  =  05h29nW9 ,6  =  -04°41'32,/7,  J2000) 
using  the  10  stations  of  the  VLBA,  which  is  operated  by  the 
National  Radio  Astronomy  Observatory  (NRAO).  Table  5 
lists  the  details  of  our  SiO  maser  observations  recorded  over 
three  epochs.  Listed  are  the  epoch,  date,  time.  Mean  Julian 
Date  (MJD)  of  the  start  of  observation,  experiment  duration, 
number  of  antennas,  and  the  visual  phase  <bvls.  The  first  two 


Table  5.  VLBA  observations  of  S  Ori. 


Ep. 

Date 

Time 

UTC 

MJD 

-2.4E6 

Dur. 

(h) 

No. 

Ant. 

®viS 

A 

2005-01-17 

02:36 

53387.6 

5.0 

9 

0.46 

B 

2005-02-28 

23:47 

53430.5 

2.6 

10 

0.56 

C 

2005-11-05 

07:23 

53679.8 

5.0 

10 

1.14 

epochs  occurred  near  stellar  minimum  and  were  spaced  42  days 
apart  or  approximately  10%  of  the  stellar  period.  The  third 
epoch  occurred  roughly  9  months  later  just  after  maximum  on 
the  next  stellar  cycle.  Data  were  recorded  in  a  similar  manner 
for  all  three  epochs  at  reference  frequencies  of  43.122080GHz 
and  42.820587  GHz  for  the  v  —  1  and  v  =  2  SiO  transitions, 
respectively.  Right-  and  left-hand  circular  polarization  was 
used  with  8-MHz  (56.1  km  s-1)  bands  centered  on  the  local 
standard-of-rest  (LSR)  velocity  of  1 8.0  kins  '.  Interspersed 
with  observations  of  S  Ori  were  10-min  scans  of  extragalactic 
continuum  sources  (0423-013  and  0359+509)  for  the  purpose 
of  delay  and  bandpass  calibration.  For  epoch  A  only  9  of  the 
10  stations  of  the  VLBA  were  available  for  our  observations. 
Unfortunately,  the  one  antenna  that  was  unavailable  was  the 
station  at  Mauna  Kea,  HI.  This  station  when  combined  with  the 
St.  Croix,  VI  antenna  provides  the  longest  baseline  (greatest 
resolution)  in  the  array.  For  this  reason  the  resolution  of  the 
epoch  A  images  is  reduced  from  that  of  the  other  two  epochs. 
There  was  also  a  problem  with  the  second  epoch  (epoch  B)  of 
observations  in  that  roughly  half  of  the  data  were  accidentally 
lost  after  the  correlation.  Despite  this  loss,  we  were  still  able 
to  use  the  remaining  data  to  produce  quality  images  of  the 
masers.  NRAO  kindly  scheduled  a  third  epoch  of  observations 
(epoch  C)  to  make  up  for  the  loss. 

4.1.  VLBA  data  reduction  and  analysis 

The  data  were  correlated  at  the  VLBA  correlator  operated  by 
NRAO  in  Socorro,  New  Mexico.  Auto  and  cross-correlation 
spectra  consisting  of  256  channels  with  channel  spacings 
of  31.25  kHz  (~0.2kms~1)  were  produced  by  the  correla¬ 
tor.  Calibration  was  performed  using  the  Astronomical  Image 
Processing  System  (AIPS)  maintained  by  NRAO.  The  data  were 
calibrated  in  a  manner  similar  to  what  was  performed  in  BW05. 
For  each  epoch,  residual  delays  due  to  the  instrumentation  were 
corrected  by  performing  a  fringe  fit  on  the  continuum  calibrator 
scans.  Residual  group  delays  for  each  antenna  were  determined 
and  applied  to  the  spectral  line  data.  Variations  in  the  residual 
delays  ranged  from  2-4  ns  resulting  in  phase  errors  of  no  more 
than  1.5-3°  across  the  8-MHz  band. 

The  bandpass  response  was  determined  from  scans  on  the 
continuum  calibrators  and  used  to  correct  the  target-source  data. 
The  time-dependent  gains  of  all  antennas  relative  to  a  reference 
antenna  were  determined  by  fitting  a  total-power  spectrum  (from 
the  reference  antenna  with  the  target  source  at  a  high  elevation) 
to  the  total  power  spectrum  of  each  antenna.  The  absolute  flux 
density  scale  was  established  by  scaling  these  gains  by  the  sys¬ 
tem  temperature  and  gain  of  the  reference  antenna.  Errors  in 
the  gain  and  pointing  of  the  reference  antenna  and  atmospheric 
opacity  variations  contribute  to  the  error  in  the  absolute  ampli¬ 
tude  calibration,  which  is  accurate  to  about  15-20%. 

Residual  fringe  rates  were  obtained  by  fringe-fitting  a 
strong  reference  feature  in  the  spectrum  of  each  maser  tran¬ 
sition.  For  epoch  A  we  used  the  channel  at  velocity  Vlsr  = 
12.3  km  s  1  for  the  u  =  I  transition  and  the  channel  at  veloc¬ 
ity  Plsr  =  16.9  km  s  1  for  the  is  =  2  transition.  For  epoch  B  we 
used  channels  at  Vlsr  =  12.8  km  s-1  and  Vlsr  =  12.1  km  s-1 


202 


M.  Wittkowski  et  al. :  The  Mira  variable  S  Orionis:  photosphere,  molecular  layer,  dust  shell,  and  SiO  maser  shell 


Fig.  7.  The  (red)  v  =  2,  J  =  1-0  (42.8GHz)  and  (green)  v  =  1,  J  =  1-0  (43.1  GHz)  maser  images  overlaid  onto  pseudo-color  representations  of 
the  infrared  intensity.  The  continuum  photosphere  (in  fact  mostly  hidden  behind  the  molecular  atmosphere  in  the  IV-band)  is  enhanced  to  a  light 
blue  color.  The  darker  blue  shades  represent  our  model  intensity  profile  as  in  panels  (e)  of  Figs.  2-5,  and  the  green  shades  represent  the  location 
of  the  AFO3  dust  shell  on  top  of  the  low-intensity  extended  wings  of  the  molecular  atmosphere.  The  true  location  of  the  star  relative  to  the  maser 
images  is  unknown.  Here,  we  assume  that  the  center  of  the  star  coincides  with  the  center  of  the  maser  spot  distribution.  Synthesized  beam  sizes  for 
epoch  A  are  0.59  x  0.40  mas  and  0.60  x  0.37  mas  for  the  v  =  1  and  v  =  2  transitions,  respectively.  For  epoch  B  these  sizes  are  0.52  x  0.21  mas  and 
0.51  x  0.20  mas,  and  for  epoch  C  0.46  x  0.16  mas  and  0.44  x  0. 16  mas,  respectively.  Epoch  D  is  a  MIDI-only  epoch  and  SiO  maser  observations 
have  not  been  obtained  at  this  epoch.  The  size  of  each  panel  is  30  x  30  mas,  corresponding  to  14.4  x  14.4  AU. 


for  the  v  —  1  and  v  =  2  transitions,  respectively.  Finally,  for 
epoch  C  we  used  the  same  velocity  channel  for  both  transi¬ 
tions  Vlsr  =  10.3  km  s  '.  The  resulting  fringe-rate  solutions 
were  applied  to  all  channels  in  each  spectrum.  An  iterative  self¬ 
calibration  and  imaging  procedure  was  then  performed  to  map 
this  reference  channel  for  each  transition.  The  resulting  residual 
phase  and  amplitude  corrections  from  the  reference  channels  at 
42.8  GHz  and  43.1  GHz  were  applied  to  all  channels  in  the  re¬ 
spective  bands.  This  procedure  was  repeated  for  each  transition 
at  each  epoch. 

In  order  to  accurately  compare  the  distributions  of  the  two 
maser  transitions,  it  is  desirable  to  determine  a  common  spa¬ 
tial  reference  point.  However,  after  the  fringe-fitting  step  to 


determine  residual  fringe-rates,  all  absolute  position  information 
is  lost  for  the  VLBA  data.  To  accomplish  this  registration,  we 
used  the  same  technique  employed  in  BW05  with  the  residual 
fringe-rates,  phases,  and  amplitudes  from  one  transition  (e.g., 
v  =  2,  42.8  GHz)  used  to  calibrate  the  other  transition  (e.g., 
v  =  1,  43.1  GHz)  in  order  to  determine  an  offset  between  the 
two  phase  centers. 

The  final  image  cubes  of  the  SiO  maser  emission  consist¬ 
ing  of  1024  x  1024  pixels  ( — 5 1  x  51  mas)  were  generated  for 
each  transition  at  each  of  the  three  epochs.  Images  were  pro¬ 
duced  for  all  spectral  channels  containing  SiO  maser  emission  in 
each  transition  at  each  epoch.  Synthesized  beam  sizes  used  for 
epoch  A  were  0.59  x  0.40  mas  and  0.60  x  0.37  mas  for  the  v  —  1 
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and  u  =  2  transitions,  respectively.  For  epoch  B  these  sizes  were 
0.52x0.21  mas  and  0.51  x  0.20  mas.  For  epoch  C  the  beam  sizes 
used  in  the  imaging  were  0.46  x  0.16  mas  and  0.44  X  0.16  mas. 
Note  that  the  resolution  was  greatly  reduced  in  epoch  A  due  to 
the  loss  of  the  Mauna  Kea,  HI  VLBA  antenna  mentioned  previ¬ 
ously.  The  resolution  of  epoch  B  is  slightly  worse  than  that  of 
epoch  C  due  to  the  loss  of  half  the  data,  thus  reducing  the  uv 
coverage  of  the  array.  Using  the  procedure  described  above,  the 
image  cubes  from  the  two  transitions  were  spatially  aligned  for 
each  epoch.  Off-source  rms  noise  estimates  in  the  images  ranged 
from  5  mJy  to  9  mJy  for  epoch  A,  3  mJy  to  7  mJy  for  epoch  B, 
and  7  mJy  to  18  mJy  for  epoch  C.  Figure  7  shows  the  total  in¬ 
tensity  images  of  the  v  =  1,  43.1-GHz  and  v  —  2,  42.8-GHz 
SiO  maser  emission  toward  S  Ori  overlaid  onto  a  pseudo-color 
representation  of  the  mid-infrared  intensity  profile  obtained  in 
Sect.  3. 

The  six  resulting  image  cubes  (2  transitions  at  3  epochs) 
were  then  analyzed  to  extract  relevant  maser  parameters.  Two- 
dimensional  Gaussian  functions  were  fit  to  maser  emission 
above  a  cut-off  flux  density  of  100  mJy  in  each  spectral  (veloc¬ 
ity)  channel  using  the  AIPS  task  SAD.  This  fitting  yielded  veloc¬ 
ities,  flux  densities,  and  positions  in  right  ascension  and  declina¬ 
tion  for  all  emission  components  identified  in  the  images.  Errors 
in  right  ascension  and  declination  of  identified  components  were 
computed  using  the  fitted  source  size  divided  by  twice  the  signal- 
to-noise  ratio  (s/n)  in  the  image  and  ranged  from  1  /.(as  for 
features  with  high  s/n,  to  28  //as  for  features  with  lower  s/n. 

The  remaining  analysis  of  the  maser  component  identifica¬ 
tions  was  performed  outside  of  the  AIPS  package.  For  compar¬ 
ison  with  the  total  flux  densities  derived  from  the  MIDI  data, 
we  summed  the  flux  densities  and  their  errors  for  all  fitted  com¬ 
ponents  in  all  channels  for  each  transition/epoch.  These  inte¬ 
grated  flux  densities  are  discussed  further  in  Sect.  5.  Since  the 
-0.2  km  s  _1  channel  spacing  is  sufficient  to  resolve  the  masers 
spectrally,  features  typically  appear  in  multiple  adjacent  spectral 
channels.  Positions  in  right  ascension  and  declination  and  center 
velocities  for  the  masers  were  determined  using  a  flux-density- 
squared  weighted  average  for  features  identified  in  two  or  more 
adjacent  channels  with  a  spatial  coincidence  of  0.5  mas.  The  flux 
assigned  to  the  maser  averages  was  the  maximum  single-channel 
flux  density.  The  results  from  the  image  analysis  and  component 
averaging  are  described  in  the  next  section. 

4.2.  VLBA  results 

4.2.1 .  The  spatial  structure  of  the  SiO  masers 
The  maser  features  characterized  using  the  procedures  described 
above  are  listed  in  Tables  A.2-A.4  and  are  shown  in  Fig.  8, 
where  there  are  six  panels  ordered  by  epoch  left  to  right  and 
by  transition  top  to  bottom.  Within  each  panel  there  are  two  sub¬ 
panels.  Each  upper  sub-panel  shows  the  spectrum  of  the  maser 
emission  with  flux  density  plotted  as  a  function  of  local  standard 
of  rest  (LSR)  velocity.  Point  sizes  are  proportional  to  the  loga¬ 
rithm  of  the  flux  density  and  are  color-coded  according  to  ve¬ 
locity  bin  in  increments  of  1 .7  km  s'1 .  The  bottom  sub-panels  of 
Fig.  8  show  the  spatial  distribution  of  the  SiO  masers  as  observed 
in  the  images.  Point  sizes  are  again  proportional  to  the  logarithm 
of  the  flux  density  with  the  same  velocity  color-coding  as  in  the 
top  sub-panels. 

Examining  the  upper  sub-panels  and  the  component  tables, 
we  find  that  the  masers  have  shifted  slightly  toward  the  blue  end 
of  the  spectrum  from  our  previous  BW05  observations.  There  the 
maser  velocities  ranged  from  1 1 . 1  to  24.7  km  s'! ,  while  for  these 


latest  three  epochs  the  masers  span  a  velocity  range  from  8.6 
to  19.0km s-1.  The  velocity  range  for  our  latest  three  epochs  is 
consistent  with  earlier  spectral  measurements  by  Jewell  et  al. 
(1991),  who  measured  a  velocity  range  of  7.1  to  19.6kms_1 
for  the  v  —  1  masers.  Spectral  ranges  measured  by  Cho  et  al. 
(1996)  were  much  narrower  at  11.0— 16. 0km  s  for  the  u  =  1 
and  12. 0—14.' 7km  s~'  for  the  v  -  2  masers,  respectively.  For  the 
three  epochs  discussed  here,  flux  densities  of  the  peak  compo¬ 
nent  in  the  spectrum  ranged  from  3.0  to  5.4  Jy  beam  1  for  the 
v  —  1  masers  and  1.3  to  6.4  Jy  beam'1  for  the  v  —  2  masers. 
These  are  close  to  our  BW05  results  for  which  we  measured 
peaks  of  4.3  Jy  beam'1  and  5.3  Jy  beam'1  for  the  v  —  1  and 
v  =  2  masers,  respectively. 

When  compared  to  the  spectrally-averaged  total-intensity 
images  (Fig.  7),  the  maser  component  maps  in  lower  sub-panels 
of  Fig.  8  provide  an  accurate  representation  of  the  spatial  struc¬ 
ture  of  the  maser  emission  summed  over  all  velocity  chan¬ 
nels  in  the  image  cube.  The  component  maps  all  show  that  the 
SiO  masers  in  both  transitions  form  partial  to  full  rings  of  emis¬ 
sion  with  the  typical  clumpy  distribution.  The  top  three  panels 
show  the  v  =  1,7  =  1-0  to  form  a  partial  ring  for  epoch  A 
with  a  gap  to  the  north,  nearly  a  full  ring  for  epoch  B,  and  again 
a  partial  ring  dominated  by  emission  in  the  northern  and  west¬ 
ern  regions  for  epoch  C.  For  the  v  =  2,7  =  1-0  SiO  masers 
shown  in  the  bottom  three  panels,  we  observe  nearly  a  full  ring 
of  emission  for  epoch  A,  a  partial  ring  with  a  gap  to  the  south 
for  epoch  B,  and  a  partial  ring  dominated  by  northern  and  west¬ 
ern  features  for  epoch  C.  In  all  three  epochs  the  v  —  1  masers 
appear  to  be  distributed  over  a  larger  area  on  the  sky  than  the 
v  —  2  masers.  In  addition,  the  v  =  1  rings  appear  thicker  than 
the  corresponding  v  =  2  rings  at  each  epoch.  It  should  be  noted 
that,  while  epochs  A  and  B  were  spaced  only  42  days  apart,  there 
are  significant  differences  in  the  number  and  distribution  of  the 
masers  in  both  transitions.  This  is  not  unexpected  since  previous 
multi-epoch  VLBI  studies  (e.g.  Diamond  &  Kemball  2003)  have 
shown  the  fine  structure  of  the  SiO  masers  to  vary  significantly 
over  periods  as  short  as  two  weeks.  This  variation,  combined 
with  the  relative  complexity  of  the  maser  structures,  precluded 
us  from  tracking  individual  maser  components  in  order  to  deter¬ 
mine  maser  proper  motions. 

To  better  characterize  the  apparent  size  and  thickness  of  the 
SiO  maser  distributions,  we  computed  the  average  angular  dis¬ 
tance  of  the  masers  from  the  center  of  their  distribution  {?).  To 
accomplish  this,  we  first  determined  the  center  of  the  distribu¬ 
tion  by  performing  a  least-square  fit  of  a  circle  to  the  combined 
v  =  1  and  v  —  2  maser  positions.  We  were  able  to  combine  the 
two  sets  of  masers  only  because  we  had  previously  aligned  the 
images  of  the  two  transitions  following  the  procedure  described 
in  Sect.  4.1.  This  fit  produced  a  common  center  from  which  we 
computed  the  mean  maser  angular  distance  ?sio  and  the  standard 
deviation  for  each  transition  independently  for  each  epoch. 

The  mean  angular  distances  from  the  center  for  the  observed 
SiO  masers  at  v  —  1,  43.1  GHz  and  v  =  2,  42.8  GHz  are  listed  in 
Table  6.  These  distances  are  also  indicated  by  a  dashed  circle  in 
each  of  the  six  panels  in  Fig.  8.  The  standard  deviations  in  the 
mean  distances  are  also  listed  in  Table  6  and  represented  in  Fig.  8 
as  the  error  on  the  mean.  Two  times  the  standard  deviation  also 
provides  an  indication  of  the  thickness  of  the  shell.  The  listed 
parameters  in  Table  6  are  for  each  epoch  of  observation/visual 
phase  and  each  of  the  two  transitions,  the  apparent  mean  angular 
distances  from  the  center  of  the  masers  spots  rv=  12,  the  ratio  of 
the  mean  distances  of  the  two  transitions  rv=i/rv=2,  as  well  as  the 
mean  LOS  velocities  Vv=it2  and  their  ratio  Vv=i/Vv=2. 
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Fig.  8.  LOS  velocity  structure  of  the  (top)  v  =  1  ,J=  1-0  (43.1GHz)  and  (bottom)  v  =  2,  J  =  1-0  (42.8  GHz)  SiO  maser  emission  toward 
S  Ori.  The  epochs  of  observations  are  from  left  to  right  Epoch  A  (Ovis  =  0.46),  Epoch  B  (<f>vis  =  0.56),  and  Epoch  C  (<Fvis  =  1.14).  The  top 
sub-panels  of  each  epoch  and  transition  show  the  spectra  formed  by  plotting  maser  intensity  versus  velocity,  color-coded  in  1.7  km  s-1  velocity 
increments  from  redward  to  blueward.  The  solid  lines  in  the  top  sub-panels  represent  the  scalar-averaged  cross-power  spectra  averaged  over  all  of 
the  VLBA  antennas.  The  bottom  (main)  sub-panels  plot  the  spatial  and  velocity  distributions  of  the  masers.  The  color  of  each  point  represents  the 
corresponding  velocity  bin  in  the  spectrum,  and  the  size  of  each  point  is  proportional  to  the  logarithm  of  the  flux  density.  Errors  in  the  positions 
of  the  features  are  smaller  than  the  data  points.  The  dashed  circles  are  based  on  the  mean  angular  distances  of  the  SiO  masers  from  the  centers 
of  the  distributions.  The  light  blue  filled  circles  in  the  centers  illustrate  for  comparison  the  angular  sizes  of  the  continuum  photosphere  derived  in 
Sect.  3.4.  The  true  location  of  the  star  relative  to  the  masers  is  unknown  and  here  assumed  to  coincide  with  the  center  of  the  maser  spot  distribution. 
We  use  also  panels  (e)  of  Figs.  2-5,  as  well  as  Fig.  7,  for  a  comparison  of  the  SiO  maser  ring  radii  to  the  extension  of  the  molecular  atmosphere 
and  the  inner  dust  shell  boundary. 


From  Table  6  we  see  that  r  ranges  from  8.6  to  10.2  mas. 
At  every  epoch  we  find  that  the  v  =  1  masers  have  a  greater 
mean  distance  from  center  than  the  v  =  2  masers.  This  finding 
is  consistent  with  BW05  in  which  we  found  rv=i  =  9.4  mas  and 
rv= 2  =  8.8  mas.  In  BW05,  however,  we  were  unable  to  make  a 
definitive  statement  about  the  relative  sizes  of  the  two  maser  dis¬ 
tributions  because  we  had  only  the  one  epoch  of  observations, 
and  at  this  epoch,  the  v  —  2  masers  were  primarily  confined  to 
a  small  region  on  the  northwest  side  of  the  shell.  In  our  more 
recent  epochs,  the  rings  of  emission  are  much  more  complete 
making  the  comparison  between  transitions  easier.  Although  the 
differences  in  r  for  epochs  A,  B,  and  C  are  less  than  the  standard 


deviation  (thickness  of  the  ring),  they  are  still  greater  than  the 
resolution  of  the  images  given  earlier  by  the  synthesized  beam 
sizes  of  0.6  x  0.4  mas  for  epoch  A  and  0.5  x  0.2  mas  for  epochs  B 
and  C.  The  fact  that  the  v  =  1  ring  is  consistently  larger  than  the 
v  =  2  ring  suggests  that  for  S  Ori  the  two  sets  of  masers  are  not 
co-spatial. 

The  relative  separation  of  the  v  —  1  and  v  -  2  masers 
for  S  Ori  is  also  consistent  with  recent  results  obtained  for  the 
SiO  masers  around  other  stars  including:  IRC  + 1001 1  (Desmurs 
et  al.  2000);  TX  Cam  (Desmurs  et  al.  2000;  Yi  et  al.  2005); 
o  Cet,  U  Ori,  and  R  Aqr  (Cotton  et  al.  2006).  In  addition, 
numerical  simulations  of  SiO  masers  (Humphreys  et  al.  1996; 
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Table  6.  Results  from  our  VLB  A  observations  of  the  v  =  1,  43.1  GHz 
and  v  =  2,  42.8  GHz  SiO  maser  transition  toward  S  Ori. 


Epoch 

Phase  <Dvis 

A 

0.46 

B 

0.56 

C 

1.14 

rv=  l  [mas] 

9.9  ±  1.4 

9.6  ±  1.1 

10.2  ±  1.2 

rB=2  [mas] 

8.6  ±  1.0 

9.1  ±  1.1 

8.9  ±  0.9 

fv=i/rv-2 

1.15 

1.05 

1.15 

V„=i  [km  s'1] 

12.6  ±2.4 

13.3  ±2.5 

13.0  ±  2.8 

V„=2  [kms-1] 

12.5  ±  2.4 

12.9  ±  3.2 

12.5  ±  2.9 

V„=i/V„=2 

1.01 

1.03 

1.04 

Gray  &  Humphreys  2000)  show  the  v  =  1 ,  7  =  1-0  masers 
to  occur  farther  from  the  star  than  the  v  —  2,  7  =  1—0  masers 
throughout  the  entire  pulsation.  The  simulations  of  Gray  & 
Humphreys  (2000)  also  show  that  the  mean  thickness  of  the 
v  -  1  shell  is  roughly  twice  that  of  the  v  -  2  shell  throughout  the 
stellar  cycle.  In  our  measurements,  we  find  the  v  =  2  shell  to  be 
thicker  in  the  BW05  observations,  the  u  —  1  shell  to  be  thicker 
in  epochs  A  and  C,  and  the  two  transitions  to  be  roughly  equiv¬ 
alent  in  epoch  B.  In  the  two  cases  for  which  the  v  -  1  shell  was 
thicker,  the  ratio  of  the  two  thicknesses  is  1.3— 1.4,  much  less 
than  the  2  determined  by  Gray  &  Humphreys.  The  angular  sizes 
of  the  maser  shells  are  reported  in  Table  7.  At  the  assumed  dis¬ 
tance  for  S  Ori  of  480  pc  ±  120  pc  (van  Belle  et  al.  2002),  1  mas 
is  roughly  equivalent  to  0.48  AU.  Thus  the  linear-scale  sizes  for 
maser  distributions  range  from  8.3  to  9.8  AU  and  the  maser  shell 
thickness  (twice  the  standard  deviation)  ranges  from  roughly  0.9 
to  1 .3  AU. 

4.2.2.  The  kinematics  of  the  SiO  masers 

In  BW05,  we  reported  that  the  SiO  masers  showed  no  coherent 
velocity  structure  indicative  of  global  expansion/infall  or  rota¬ 
tion.  Examining  Fig.  8,  we  again  find  no  evidence  of  such  co¬ 
herence  in  the  line-of-sight  (LOS)  velocities  of  the  masers  as  a 
function  of  spatial  location.  We  did,  however,  notice  that  there 
appears  to  be  a  velocity  gradient  at  all  epochs,  with  masers  to¬ 
ward  the  blue-  and  red-shifted  ends  of  the  spectrum  lying  closer 
to  the  center  of  the  distribution  than  masers  at  intermediate  ve¬ 
locities.  This  phenomenon  was  also  found  in  BW05  and  has  been 
previously  observed  in  other  SiO  maser  studies  (e.g.  Boboltz  & 
Marvel  2000;  Hollis  et  al.  2001). 

The  LSR  stellar  velocity  of  S  Ori  has  not  been  measured 
particularly  well.  In  the  General  Catalog  of  Stellar  Radial 
Velocities  (Wilson  1953),  the  value  listed  for  S  Ori  is  22kms-1. 
Young  (1995)  observed  the  sub-millimeter  CO  (3-2)  and 
CO  (4-3)  lines  toward  S  Ori  and  determined  values  of  Vlsr  = 
14.5  +  0.2  and  14.1  +  0.5  km  s-1  from  fits  to  the  spectra,  re¬ 
spectively.  Winters  et  al.  (2003)  derived  Vlsr  =  14.0km  s-1 
based  on  fits  of  the  CO  (1-0)  and  CO  (2-1)  lines.  Since  the 
SiO  masers  are  not  distributed  evenly  about  any  of  the  above  ve¬ 
locities,  we  computed  an  average  LOS  velocity  from  the  masers 
themselves  for  each  transition  at  epoch.  The  values  of  VLOS  for 
each  of  the  epochs  are  reported  in  Table  6  with  errors  of  15% 
estimated  from  the  uncertainty  of  the  absolute  amplitude  cali¬ 
bration.  Although  the  differences  are  well  within  the  error  bars, 
it  appears  that  the  v  -  2  masers  are  shifted  slightly  toward 
the  blue  end  of  the  spectrum  relative  to  the  v  =  1  masers  at 
all  three  epochs.  In  addition  to  computing  values  for  individual 
transitions  and  epochs,  we  also  combined  the  three  epochs  and 
computed  a  mean  LOS  velocity  for  each  transition.  We  obtained 
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Fig.  9.  Histogram  of  the  LOS  velocities  of  our  SiO  maser  spots  com¬ 
bining  VLB  A  epochs  A-C.  Panel  a)  shows  the  histogram  for  the  v  =  2 
42.8  GHz  masers  and  panel  b)  for  the  v  =  1  43.1  GHz  masers.  Also  indi¬ 
cated  are  the  mean  values  V„=i  and  V„=2.  K=i  is  slightly  redder  than  V„=2. 


mean  values  of  Vlsr  =  12.6  +  2.8kms~'  and  13.0  +  2.6  km  s-1 
for  the  42.8  GHz  and  43.1  GHZ  masers,  respectively.  These  val¬ 
ues  are  shown  in  Fig.  9.  Again  we  see  that  V„=1  is  slightly  redder 
than  Vd=2. 

To  verify  the  notion  that  the  higher- velocity  masers  are  found 
closer  to  the  star,  we  plotted  component  distance  from  the  center 
of  the  distribution  (radius)  versus  its  velocity  relative  to  the  av¬ 
erage  LOS  velocity  for  that  transition  (Vlos  -  Vlos)-  The  results 
are  shown  in  Figs.  10  and  1 1  for  the  v  =  1,7=  1-0  43.1  GHz 
and  the  v  -  2,  J  -  1-0  42.8  GHz  transitions,  respectively.  In 
the  figures  we  plotted  each  epoch  with  a  different  symbol.  The 
masers  in  both  transitions  appear  to  have  the  same  distribution 
with  a  central  peak  near  Vlos  -  Vlos  =  0  and  decreasing  maser 
radius  as  a  function  of  increasing  values  of  |  Vlos  -  VlosI- 

A  simple  model  that  is  often  used  in  the  case  of  OH  (e.g. 
Reid  et  al.  1977)  and  H20  (e.g.  Yates  &  Cohen  1994)  maser 
kinematics  is  that  of  a  uniformly  expanding  thin  shell.  In  this 
model  the  projected  distance  from  center  (a)  of  a  maser  on  the 
shell  is  related  to  its  LOS  velocity  V  -  V*  by  the  expression: 


(V-V+ 
\  Vsh 


=  1, 


(6) 


where  rs|,  and  Vsh  are  the  radius  and  the  expansion  velocity  of  the 
shell,  respectively.  The  LSR  velocity  of  the  star  (V*)  is  assumed 
here  to  be  equal  to  Vlos-  Numerical  simulations  of  SiO  masers 
often  employ  the  Sobolev  or  large  velocity  gradient  (LVG)  ap¬ 
proximation  (Lockett  &  Elitzur  1992;  Bujarrabal  1994;  Doel 
et  al.  1995;  Humphreys  et  al  2002).  The  LVG  approximation  as¬ 
sumes  a  large  velocity  gradient  across  the  SiO  maser  region,  thus 
producing  large  Doppler  shifts  that  serve  to  disconnect  a  region 
of  velocity  coherence  from  the  surrounding  medium  (Bujarrabal 
1994).  The  parameter  typically  used  to  characterize  the  velocity 
field  in  the  maser  region  is  the  logarithmic  velocity  gradient  e 
given  by 

d  In  V  rdV 

e  =  — —  =  -  —  •  (7) 

dlnr  V  dr 
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Fig.  10.  Distances  of  the  43.1  GHz  maser  components  from  the  com¬ 
mon  center  of  the  distribution  (radius)  versus  their  velocity  relative  to 
the  average  LOS  velocity  for  the  43.1GHz  transition  (Vlos  -  Vlos)- 
Maser  components  of  the  different  epochs  are  denoted  by  different  sym¬ 
bols  and  colors.  The  plot  indicates  that  higher- velocity  masers  are  found 
closer  to  the  star.  This  trend  can  be  explained  by  means  of  expanding 
maser  shells.  The  curves  represent  shells  at  the  inner  and  outer  bound¬ 
aries  of  the  SiO  maser  shell  using  (solid  lines)  a  logarithmic  velocity 
gradient  and  (dashed  lines)  constant  velocity  expansion.  See  text  for 
details. 


A  value  of  e  —  0  is  equivalent  to  constant  velocity  expansion, 
while  e  >  1  corresponds  to  a  velocity  field  with  large  radial  ac¬ 
celerations.  The  SiO  maser  numerical  simulations  (e.g.  Lockett 
&  Elitzur  1992;  Bujarrabal  1994;  Doel  et  al.  1995)  typically  use 
a  value  of  e  =  1.  Chapman  &  Cohen  (1986)  determined  values 
of  e  for  the  various  maser  species  surrounding  the  supergiant 
VX  Sgr  and  found  e  *  1  in  the  SiO  region,  e  ~  0.5  in  the  region 
of  the  H2O  and  mainline  (1665  MHz,  1667  MHz)  OH  masers, 
and  e  ~  0.2  in  the  1612  MHz  OH  maser  region. 

Plotted  along  with  the  data  in  Figs.  10  and  1 1  are  four  curves 
computed  using  the  above  equations.  They  represent  shells  at  the 
inner  and  outer  boundaries  of  the  SiO  masers.  For  the  masers 
in  each  transition,  we  used  an  inner  shell  radius  equivalent  to 
the  smallest  average  distance  minus  its  standard  deviation  from 
Table  6.  Similarly,  for  the  outer  shell  radius  we  used  the  largest 
average  distance  plus  its  standard  deviation  from  Table  6.  For 
the  scenario  of  a  logarithmic  velocity  gradient,  the  velocity  at 
the  inner  v  —  2  boundary  was  set  at  7kms-1  and  the  velocities 
of  the  other  shells  computed  assuming  e  =  1 .  The  velocity  of  the 
outermost  shell  of  the  v  =  1  transition  is  Lh  ~  10.5  km  s  1 . 

For  the  scenarion  of  a  constant  velocity  expansion  (e  =  0), 
the  best  by-eye  fit  results  in  Vsh  =  lOkms-1.  In  this  case,  the 
higher  expansion  velocity  at  the  inner  boundary  of  the  v  —  2 
transition  is  required  in  order  to  provide  a  reasonable  fit  to  the 
outer  boundaries  of  both  transitions.  For  the  e  —  1  case,  the  inner 
v  =  2  velocity  can  be  lower  with  the  radial  acceleration  allow¬ 
ing  for  a  better  fit  to  the  outer  boundaries.  For  comparison,  the 
escape  velocities  computed  for  the  innermost  v  =  2  and  outer¬ 
most  v  =  1  shells,  assuming  a  mass  for  S  Ori  of  1.0  Me,  are  22.1 
and  1 8.0  km  s'1,  respectively.  Therefore,  even  at  the  velocities 


Fig.  11.  As  Fig.  10,  but  for  the  42.8  GHz  transition.  The  axes’  scales  are 
identical  for  the  two  figures.  Compared  to  Fig.  10,  it  is  evident  that  the 
42.8  GHz  maser  spots  lie  at  systematically  closer  distances  to  the  center 
of  the  common  maser  distribution  than  the  43. 1  GHz  maser  spots. 


bracketed  by  the  curves,  the  maser  gas  is  still  gravitationally 
bound. 

From  the  SiO  maser  data  on  S  Ori,  we  cannot  definitively 
rule  out  constant  velocity  expansion;  however,  the  SiO  maser  re¬ 
gion  is  very  dynamic,  influenced  by  the  pulsation  of  the  star  and 
the  passage  of  shocks.  Theoretical  studies  all  seem  to  include 
some  form  of  velocity  gradient  across  the  SiO  maser  region.  In 
addition,  constant  velocity  expansion  would  be  more  conducive 
to  masers  along  the  line  of  sight  to  the  star  because  of  the  in¬ 
creased  velocity  coherence.  SiO  masers  are  rarely  observed  in 
front  of  the  stellar  disc,  and  for  S  Ori  in  particular,  such  masers 
have  not  be  seen.  We  therefore  consider  some  form  of  accelera¬ 
tion  in  the  region  the  more  likely  scenario. 

The  canonical  interpretation  of  the  rings  formed  by  the 
SiO  masers  is  that  the  masers  are  tangentially  amplified  and  pri¬ 
marily  confined  to  a  region  close  to  the  plane  of  the  sky.  Under 
this  interpretation,  the  ring  size  and  thickness  are  direct  indica¬ 
tors  of  the  equivalent  3-D  shell  size  and  thickness.  The  radially 
expanding  spherical  shell  described  above  will  naturally  result  in 
a  projection  effect  where  masers  out  of  the  plane  of  the  sky  have 
smaller  angular  distances  from  center  due  to  the  projection.  In 
the  case  of  S  Ori,  this  projection  is  rather  small.  The  equation  for 
the  angle  8  measured  between  the  LOS  and  a  point  on  the  shell 
locating  the  maser  can  be  written  in  terms  of  the  velocities  as 


sin  8  = 


(8) 


In  the  most  extreme  case,  all  of  the  thickness  observed  in  the  ring 
is  due  to  projection  of  the  masers  on  the  outermost  thin  shell  with 
a  velocity  of  Ush  ~  1 0.5  km  s~'  (e  =  1).  Figures  9-11  show  that 
the  majority  of  the  masers  are  within  ~4  k m  s  1  of  LOS  cen¬ 
ter  velocity.  For  these  extreme  4km s_1  masers,  sin  8  =3  0.92 
and  8  ~  67°.  This  results  in  a  difference  of  ~8%  between  the 
projected  distance  of  the  maser  and  its  true  radial  distance.  For 
masers  less  than  4  km  s  1  the  difference  is  smaller.  Assuming 
the  masers  are  distributed  uniformly,  the  differences  between  the 
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Table  7.  Overview  of  the  photometric  and  spatial  parameters  of  S  Ori 
as  a  function  of  stellar  phase,  derived  from  both  our  VLTI  and 
VLB  A  observations. 


A 

B 

C 

D 

Phase  4>vis 

0.44 

0.55 

1.15 

1.27 

V  [mag] 

12.8  ±0.3 

14.0  ±  0.3 

9.4  ±  0.3 

10.5  ±  0.4 

fN  [Jy] 

187  ±  12 

152  ±  17 

198  ±  18 

203  ±  13 

/43.1  [Jy] 

119  ±  18 

203  ±  30 

198  ±  30 

/42.8  [Jy] 

148  ±  22 

46  ±7 

383  ±  57 

©Phot  [mas] 

9.0  ±  0.3 

7.9  ±0.1 

9.7  ±0.1 

9.5  ±  0.4 

©in  [mas] 

16.2  ±  1.9 

15.9  ±  1.6 

21.2  ±  1.9 

22.8  ±  2.1 

©4i  t  Tmas] 

19.6  ±  2.8 

19.0  ±  2.6 

20.4  ±  2.4 

©42.8  [mas] 

19.0  ±  2.0 

17.8  ±  3.4 

17.8  ±  1.8 

©in /©Phot 

1.8  ±  0.2 

2.0  ±  0.2 

2.2  ±  0.2 

2.4  ±  0.2 

©43.1 /©Phot 

2.2  ±  0.3 

2.4  ±  0.3 

2.1  ±0.3 

©42.8 /©Phot 

2.1  ±0.2 

2.3  ±  0.4 

1.9  ±0.2 

©43.1 /©in 

1.2  ±0.2 

1.2  ±  0.2 

1.0  ±0.1 

©42.8 /©in 

1.2  ±0.2 

1.1  ±  0.2 

0.8  ±0.1 

J^Phot  [Re] 

469  ±  127 

411  ±  111 

498  ±  134 

491  ±  133 

projected  ring  sizes  and  thicknesses  and  their  3-D  equivalents 
is  -4%.  Thus,  the  given  errors  of  the  projected  ring  sizes  deter¬ 
mined  from  the  maser  distributions  include  the  uncertainty  of  the 
projection  effect. 

5.  VLTI/MIDI  and  VLBA  results  as  a  function 
of  stellar  phase 

Table  7  provides  an  overview  on  the  photometric  and  spatial  pa¬ 
rameters  of  S  Ori  as  a  function  of  observational  epoch/stellar 
phase,  as  derived  from  both  our  VLTI/MIDI  and  VLBA  obser¬ 
vations  in  Sects.  3  and  4.  The  listed  photometric  parameters  in¬ 
clude  the  V  magnitude  based  on  the  lightcurve  in  Fig.  1,  the 
/V-band  flux  based  on  our  photometric  MIDI  data  (Sect.  3.1), 
and  the  total  flux  values  of  the  two  maser  transitions  (Sect.  4.1). 
The  listed  spatial  information  includes  the  photospheric  angu¬ 
lar  diameter  ©phot  (Sect.  3.4),  the  inner  dust  shell  diameter  0jn 
(Sect.  3.4),  and  the  maser  ring  angular  diameters  ©43.1  and  @42.8 
(Sect.  4.2.1).  To  compare  the  continuum  photospheric  diameter, 
the  inner  dust  shell  diameter,  and  the  maser  shell  diameters,  the 
ratios  ©in/©Phot,  ©43.i/42.8/©Phot,  and  ©43.i/42.8/©in  are  listed  as 
well.  As  an  absolute  scale,  the  stellar  photospheric  radius  Rpu(„ 
at  each  phase  is  derived  from  ©phot  and  the  adopted  distance  to 
S  Ori  (Sect.  2).  Here,  the  variability  phases  are  the  mean  values 
of  the  respective  phases  of  the  MIDI  and  VLBA  observations. 

The  continuum  photospheric  radii  /?ph0t,  the  inner  dust  shell 
radii  R[n  (corresponding  to  0jn)  and  the  maser  shell  radii  R43.] 
and  R42.8  (corresponding  to  ©43.1/42.8)  are  indicated  in  panels  (d) 
of  Figs.  2-5,  which  show  the  synthetic  mid-infrared  CLVs.  This 
allows  a  comparison  of  these  radii  to  each  other  and  to  the  profile 
of  the  extended  molecular  atmosphere.  Figure  7  provides  as  well 
a  graphical  comparison  of  the  photospheric  disc,  the  molecular 
layers,  the  dust  shell,  and  the  maser  spots  for  each  epoch. 

Variability  of  visual,  mid-infrared,  and  maser  flux  densities 
The  visual,  mid-infrared,  and  maser  fluxes  all  show  clear  vari¬ 
ability  as  a  function  of  observational  epoch.  The  stellar  phases 
given  throughout  this  paper  are  visual  phases  as  introduced  in 
Sect.  2.  The  V-band  flux  fy  shows  clear  variability  with  an  am¬ 
plitude  of  about  50  Jy  (-30%).  fy  exhibits  a  minimum  at  our 


epoch  B  at  visual  phase  0.55,  i.e.  shortly  beyond  the  visual 
minimum,  rises  toward  stellar  maximum,  and  continues  to  in¬ 
crease  between  epochs  C  and  D  at  visual  phases  1.15  and  1.27, 
probably  indicating  a  phase  lag  and  different  shape  of  the  N- 
band  variability  curve  compared  to  the  visual  lightcurve.  The 
total  maser  fluxes  are  significantly  lower  near  stellar  minimum 
(epochs  A  and  B)  compared  to  stellar  maximum  (epoch  C),  with 
the  exception  of  the  strong  .1  value  at  epoch  B.  The  latter  can 
most  likely  be  explained  by  the  larger  number  of  maser  spots 
forming  a  nearly  complete  ring-like  structure  for  this  transition 
and  epoch,  while  the  other  transitions  and  epochs  show  only  par¬ 
tial  maser  rings. 

Photospheric  and  dust  shell  radii  as  a  function  of  visual 
phase  The  photospheric  angular  diameter  ©phot  derived  from 
our  MIDI  observations  and  modeling  shows  a  clearly  phase- 
dependent  size  with  an  amplitude  of  about  1.8  mas  (—20%)  that 
is  well-correlated  in  phase  with  the  visual  lightcurve.  It  shows 
a  minimum  at  epoch  B,  which  is  nearest  to  the  visual  mini¬ 
mum,  and  a  maximum  at  epoch  C,  which  is  nearest  to  the  vi¬ 
sual  maximum.  Also  the  inner  dust-shell  radius  shows  a  varia¬ 
tion  (amplitude  -7  mas  -30%),  which  is  more  closely  correlated 
in  phase  and  amplitude  to  fy  than  to  V.  The  inner  dust  shell  ra¬ 
dius  is  smallest  at  our  epoch  B  near  visual  minimum  and  largest 
at  epoch  D  beyond  visual  maximum  (phase  1.27).  Compared 
to  ©phot,  the  inner  dust  shell  boundary  is  located  significantly 
closer  to  the  stellar  surface  at  our  epochs  A  and  B  near  visual 
minimum  (1. 8-2.0  x  ©phot)  than  at  the  post-maximum  visual 
phases  C  and  D  (2. 2-2.4  x  ©phot)-  Simultaneously,  as  discussed 
in  Sect.  3.4,  the  density  gradient  p  is  significantly  steeper  and 
the  optical  depth  significantly  higher  near  visual  minimum  than 
at  the  post-maximum  visual  phases.  This  can  be  understood  by 
a  wind  model  with  increased  mass  loss  near  minimum  phase  oc¬ 
curring  close  to  the  stellar  surface  and  an  expanded  dust  shell 
after  maximum  (cf.  Sect.  3.4). 

Maser  ring  radii  as  a  function  of  visual  phase  The  mean 
SiO  maser  ring  radii  obtained  at  our  three  VLBA  epochs  show  a 
small  variability  with  amplitude  of  <7%,  which  is  only  a  fraction 
of  the  mean  width  of  the  maser  shells  of  -13%  and  clearly  below 
the  variability  amplitude  of  the  photospheric  and  dust  shell  radii 
at  the  same  epochs.  However,  we  cannot  rule  out  that  the  maser 
ring  radii  show  a  greater  variability  in  between  our  epochs  of 
observation.  Observations  by  Cotton  et  al.  (2006)  consistently 
find  variations  in  maser  ring  diameters  of  3-14%  for  different 
Mira  stars,  but  no  clear  correlation  with  visual  variability  phase. 
Theoretical  estimates  by  Humphreys  et  al.  (2002)  predict  a  vari¬ 
ability  of  the  43.1GHz  maser  shell  radius  with  an  amplitude 
of  -20%. 

Maser  ring  radii  compared  to  the  photospheric  radius  The 
true  location  of  the  star  relative  to  the  maser  spots  is  unknown. 
Based  on  the  assumption  that  the  center  of  the  maser  spot  distri¬ 
bution  coincides  with  the  center  of  the  star,  a  comparison  of  the 
maser  ring  radii  with  the  continuum  photospheric  radii  gives  an 
estimate  of  the  distance  of  the  maser  spots  from  the  stellar  pho¬ 
tosphere.  The  ratios  between  the  maser  ring  radii  and  the  photo¬ 
spheric  radii  in  Table  7  range  between  2.1  +  0.3  and  2.4  +  0.3. 
These  values  are  consistent  with  the  values  reported  in  BW05 
(1.9  and  2.0),  with  those  by  Cotton  et  al.  (2004),  which  range  for 
different  Mira  stars  -  not  always  compared  at  the  same  phase  - 
between  1.8  and  2.9,  and  the  value  by  Fedele  et  al.  (2005)  for 
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R  Leo  of  1.8.  Our  observations  indicate  a  variability  in  the  ratio 
of  the  mean  maser  shell  radius  to  the  photospheric  radius  with 
an  amplitude  of  ~10%  for  the  43.1GHz  transition,  and  ~20% 
for  the  42.8  GHz  transition,  with  uncertainties  on  the  same  or¬ 
der.  This  observed  variability  is  dominated  by  the  variability  of 
the  stellar  continuum  photospheric  diameter. 

Maser  ring  radii  compared  to  the  molecular  atmosphere 

Panels  (d)  of  Figs.  2-5  as  well  as  Fig.  7  show  that  the  mean 
maser  ring  radii  at  all  our  epochs  mark  the  region  just  beyond  the 
steepest  decrease  in  the  mid-infrared  model  intensity,  i.e.  shortly 
outward  of  the  densest  region  of  the  molecular  atmosphere. 

Maser  ring  radii  compared  to  the  inner  dust  shell  bound¬ 
ary  Table  7,  as  well  as  panels  (d)  of  Figs.  2-5,  indicate  that 
the  inner  boundary  of  the  AI2O3  dust  grains  near  visual  min¬ 
imum  (epochs  A  and  B)  appears  to  be  located  between  the 
steepest  decrease  in  the  mid-infrared  model  intensity  and  the 
mean  location  of  the  SiO  maser  spots.  This  means  that  the  low- 
density  atmospheric  molecular  layers,  the  SiO  maser  spots,  and 
the  ALO3  dust  grains  are  co-located  near  visual  minimum.  At 
our  post-maximum  epochs  C  and  D,  the  inner  boundary  of  the 
ALO3  dust  shell  appears  significantly  expanded,  while  the  maser 
shells  remain  at  about  the  same  location.  As  a  result,  the  inner 
boundary  of  the  ALO?  dust  grains  appears  to  have  expanded  out¬ 
ward  of  the  mean  SiO  maser  rings  at  the  post-maximum  phase. 

6.  Summary,  conclusions,  and  discussion 

We  have  obtained  the  first  time  series  of  observations  of  a 
Mira  variable  star  that  include  concurrent  infrared  and  radio  in¬ 
terferometry.  For  the  present  paper,  we  obtained  mid-infrared 
interferometry  of  S  Ori  with  VLTI/MIDI  at  four  epochs  in 
December  2004,  February/March  2005,  November  2005,  and 
December  2005.  We  observed  v  =  1,7  =  1-0  (43.1GHz)  and 
v  =  2,7  =  1-0  (42.8  GHz)  SiO  maser  emission  toward  S  Ori 
with  the  VLBA  in  January,  February,  and  November  2005,  con¬ 
current  to  within  5%  of  the  variability  period  to  the  first  three 
VLTI/MIDI  epochs.  The  first  two  common  VLTI/VLBA  epochs 
are  located  near  visual  minimum  (phases  0.44  and  0.55)  and  the 
third  after  next  visual  maximum  (phase  1.15).  The  fourth  MIDI 
epoch  occurred  at  a  later  post-maximum  phase  of  1 .27. 

The  mid-infrared  MIDI  data  are  sensitive  to  the  structure 
of  the  atmosphere  consisting  of  the  continuum  photosphere  and 
overlying  molecular  layers,  as  well  as  to  the  properties  of  the 
dust  shell.  The  MIDI  visibilities  and  flux  densities  are  modeled 
using  the  recent  M  model  atmosphere  series  by  Ireland  et  al. 
(2004b),  which  is  a  dust-free  self-excited  dynamic  model  at¬ 
mosphere  series,  to  which  we  add  an  ad-hoc  radiative  transfer 
calculation  of  the  dust  shell.  The  constrained  model  parameters 
include  the  continuum  photospheric  angular  diameter,  as  well 
as  characteristics  of  the  dust  shell.  The  MIDI  visibility  and  to¬ 
tal  flux  data  at  all  epochs  can  be  described  well  by  the  chosen 
approach  of  combining  the  dust-free  M  model  atmosphere  se¬ 
ries  with  an  ad-hoc  radiative  transfer  model  of  the  dust  shell. 
Different  model  parameters  can  be  constrained  by  different  sub¬ 
sets  of  our  MIDI  data. 

The  resulting  continuum  photospheric  angular  diameters 
at  visual  variability  phases  0.42,  0.55,  1.16,  and  1.27  are 
9.0  +  0.3  mas,  7.9 +  0.1  mas,  9.7  ±  0.1  mas,  and  9.5  ±0.4 mas,  re¬ 
spectively.  The  photospheric  radius  thus  shows  a  significantly 
phase-dependent  size  with  amplitude  of  ~20%  that  is  well  cor¬ 
related  in  phase  with  the  visual  lightcurve. 


The  dust  shell  can  be  modeled  well  with  ALO3  grains 
alone,  and  our  data  show  no  indication  of  silicates,  con¬ 
sistent  with  the  modeling  of  the  IRAS  LRS  spectra  by 
Lorenz -Martins  &  Pompeia  (2000).  The  inner  dust  shell  angular 
diameters  are  16.2  ±  1.9  mas,  15.9  ±  1.6  mas,  21.2  ±  1.9  mas,  and 
22.8  +  2.1  mas,  respectively,  for  the  visual  phases  given  above, 
corresponding  to  1.8,  2.0,  2.2,  and  2.4  times  the  continuum  pho¬ 
tospheric  diameters.  Thus,  dust  starts  to  form  very  close  to  the 
stellar  surface,  corresponding  to  the  region  of  the  low-intensity 
wings  of  the  atmospheric  molecular  layers  of  our  dust-free  atmo¬ 
sphere  model.  In  particular,  the  inner  dust  shell  radius  is  smallest 
at  the  near-minimum  phases,  located  directly  outside  of  the  steep 
decrease  in  the  dust-free,  mid-infrared  model  intensity.  The  inner 
dust-shell  boundary  expands  to  larger  distances  from  the  stellar 
surface  at  our  post-maximum  phases.  The  overall  ALO3  opti¬ 
cal  depth,  as  well  as  the  density  gradient,  is  largest  at  our  near¬ 
minimum  phases  and  significantly  lower  at  the  post-maximum 
stellar  phases.  These  dust  shell  characteristics  can  be  explained 
by  a  radiatively  driven  wind  with  increased  mass-loss  near  vi¬ 
sual  stellar  minimum  occurring  close  to  the  stellar  surface,  and 
an  expanded  dust  shell  at  the  post-maximum  phases. 

The  43.1GHz  and  42.8  GHz  maser  spots  show  the  typi¬ 
cal  structure  of  partial  to  full  rings  with  a  clumpy  distribu¬ 
tion.  There  is  no  sign  of  a  globally  asymmetric  gas  distribu¬ 
tion.  The  43.1  GHz  mean  maser  ring  angular  diameters  at  visual 
phases  0.46,  0.56,  and  1.14  are  19.6  ±2.8  mas,  19.0  ±2.6  mas, 
and  20.4  ±2.4  mas,  respectively.  The  values  for  the  42.8  GHz 
maser  are  systematically  smaller  by  3-19%  and  amount  to 
19.0  ±2.0  mas,  17.8  ±3. 4 mas,  and  17.8  ±1.8  mas.  The  indica¬ 
tion  of  systematically  smaller  mean  maser  ring  diameters  of 
the  42.8  GHz  transition  compared  to  the  43.1  GHz  transition  is 
consistent  with  earlier  results  in  the  literature  on  other  Mira 
stars.  The  velocity  structure  of  both  maser  transitions  indicate 
an  expansion  of  the  maser  shell.  This  expansion  is  most  likely 
an  accelerated  expansion,  but  we  cannot  completely  rule  out  a 
uniform  expansion.  The  expansion  velocity  for  both  transitions 
ranges  between  about  7  km  s-1  and  about  10  km  s_1  between  the 
inner  and  outer  shell  radii  based  on  a  scenario  of  an  accelerated 
expansion  or  amounts  to  about  10km s_1  based  on  a  scenario 
with  uniform  expansion. 

Compared  to  the  continuum  photospheric  radius  Rphot- 
the  43.1  GHz/42.8  GHz  SiO  maser  ring  radii  are  located  at 
2.2  ±0.3/2. 1+0.2  (phase  0.44),  2.4  +0.3/2.3  ±0.4  (0.55),  and 
2.1  ±0.3/1. 9  +  0.2  (1.15)  times  Rphot-  These  values  are  consis¬ 
tent  with  theoretical  estimates  and  earlier  observations  available 
in  the  literature.  Compared  to  the  molecular  atmosphere,  the 
maser  shells  mark  the  region  of  the  molecular  atmosphere  just 
beyond  the  steepest  decrease  in  the  mid-infrared  model  intensity. 
Compared  to  the  dust  shell,  the  maser  spots  are  co-located  with 
the  inner  dust  shell  near  stellar  minimum.  At  our  post-maximum 
phases,  the  maser  spots  remain  at  about  the  same  location,  while 
the  dust  shell  has  expanded  outward. 

Altogether,  our  observations  indicate  a  pulsation  of  the 
stellar  photosphere  approximately  in  phase  with  the  visual 
lightcurve.  Our  measurements  suggest  an  increased  mass-loss 
rate  near  stellar  minimum  and  formation  of  ALO3  dust  grains 
occurring  close  to  the  stellar  surface  at  about  1. 8-2.0  pho¬ 
tospheric  radii  directly  outward  of  the  dense  molecular  at¬ 
mosphere,  co-located  with  the  extended  low-intensity  wings 
of  the  molecular  atmosphere  and  with  the  SiO  maser  spots. 
At  the  post-maximum  phases,  the  dust  shell  has  expanded  to 
about  2. 2-2.4  photospheric  radii,  while  the  mean  radius  of  the 
maser  spots  remains  at  ~2.0  photospheric  radii. 
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Fig.  12.  Sketch  of  the  radial  structure  of  S  Ori's  CSE  at  {left)  near-minimum  and  (right)  post-maximum  visual  phase  as  derived  in  this  work.  Shown 
are  the  locations  of  the  continuum  photosphere  (dark  gray),  the  at  A-band  optically  thick  molecular  atmosphere  (medium  dark  gray),  the  at  (V-band 
optically  thin  molecular  atmosphere  (light  gray),  the  AI2O3  dust  shell  (dashed  arcs),  and  the  42.8  GHz  and  43. 1  GHz  maser  spots  (circles/triangles). 
The  numbers  below  and  beside  the  panels  are  the  mean  values  of  (left)  epochs  A  &  B  and  (right)  epochs  C  &  D  from  Table  7. 


Figure  12  shows  a  sketch  of  S  Ori  and  its  CSE  near  vi¬ 
sual  minimum  and  post-maximum  illustrating  the  findings  of  the 
present  work.  Our  results  on  S  Ori  are  largely  consistent  with 
the  canonical  scheme  of  the  CSE  of  AGB  stars  discussed  for 
instance  in  Reid  &  Menten  (1997).  However,  our  results  indi¬ 
cate  an  inner  dust  shell  boundary  that  is  located  closer  to  the 
optical  photosphere  compared  to  the  Reid  &  Menten  scheme  at 
about  two  photospheric  radii,  and  SiO  maser  shells  that  are  co¬ 
located  with  the  low-intensity  molecular  atmosphere  and  with 
AFO;  dust  grains  near  visual  minimum.  In  the  case  of  S  Ori, 
the  dust  shell  consists  only  of  ALO3  grains.  For  Mira  stars  that 
also  show  a  dust  shell  of  silicates,  the  inner  boundary  of  silicate 
grains  is  still  expected  at  greater  distances  of  A4  photospheric 
radii  (Lorenz-Martins  &  Pompeia  2000).  In  addition  to  a  static 
picture  of  a  Mira  star’s  CSE,  our  observations  suggest  the  forma¬ 
tion  of  a  dust  shell  close  to  the  surface  near  minimum  phase  and 
an  expansion  by  ~20%  of  the  inner  dust-shell  boundary  toward 
post-maximum  phase. 

Our  modeling  approach  of  adding  the  CLVs  of  a  dust-free 
M  model  to  that  of  a  radiative  transfer  model  of  the  pure  dust 
shell  includes  two  uncertainties.  Firstly,  the  M  model  series  is 
designed  for  the  prototype  Miras  o  Cet  and  R  Leo  and  used  here 
for  the  longer-period  Mira  S  Ori.  Secondly,  we  omit  a  possible 
interaction  of  gas  and  dust  outward  of  the  inner  dust  shell  bound¬ 
ary.  Both  effects  might  add  some  additional  uncertainty  to  the  ab¬ 
solute  values  of  the  obtained  photospheric  diameter  values  and 
dust  shell  characteristics.  However,  our  data  can  be  described 
well  by  our  modeling  approach  and  these  additional  uncertain¬ 
ties  are  not  expected  to  be  significantly  larger  than  the  diameter 
errors  given  in  Table  7  and  the  step  sizes  of  the  dust  shell  parame¬ 
ters  in  Table  4.  Furthermore,  the  relative  accuracy  of  the  obtained 
parameters  among  the  four  epochs  of  observations  are  clearly 


less  affected  by  possible  systematic  uncertainties.  Currently,  the 
chosen  approach  appears  to  be  the  best  available  one.  Similar 
studies  of  S  Ori  and  other  long-period  Mira  stars  in  the  future 
would  benefit  from  new  dynamic  model  atmosphere  series  tai¬ 
lored  to  Mira  stars  with  a  longer  period  and  a  higher  mass  com¬ 
pared  to  o  Cet  and  R  Leo,  as  well  as  from  self-consistent  calcu¬ 
lations  of  dynamic  atmospheres  and  dust  formation. 

Our  mid-infrared  interferometric  data,  in  particular  at  the 
lower  edge  of  MIDI’s  wavelength  range,  together  with  the  use  of 
the  M  model  series  provide  constraints  on  the  continuum  pho¬ 
tospheric  angular  diameter  of  S  Ori.  The  resulting  values  are 
in  good  agreement  with  several  earlier  diameter  measurements 
obtained  in  the  near-infrared  /('-band  by  different  instruments 
and  at  different  stellar  phases,  which  increases  the  confidence  in 
our  modeling  approach.  Additional  confidence  in  the  continuum 
photospheric  diameters  and  the  validity  of  the  dynamic  model 
atmospheres  can  in  future  be  reached  by  time  series  of  high  spec¬ 
tral  resolution  (R  =  1000-10000)  near-infrared  interferometry 
with  the  VLTI/AMBER  (Petrov  et  al.  2003)  instrument  concur¬ 
rent  to  VLTI/MIDI  observations. 

Our  results  obtained  at  4  epochs  within  one  pulsation  cycle 
suggest  the  formation  of  a  dust  layer  consisting  of  ALO3  grains 
close  to  the  stellar  surface  near  stellar  minimum  and  an  expanded 
dust  shell  after  stellar  maximum.  Our  understanding  of  the  con¬ 
nection  between  the  dust  formation  frequency  and  the  stellar  pul¬ 
sation  will  be  further  improved  by  using  more  epochs  within  one 
pulsation  cycle.  Dust-forming  episodes  may  also  have  their  own 
timescale  that  may  not  be  directly  related  to  the  visual  variabil¬ 
ity  (e.g.,  Hofner  &  Dorfi  1997;  Nowotny  et  al.  2005),  a  study  of 
which  requires  observations  extended  over  more  than  one  pul¬ 
sation  cycle.  Since  the  data  of  S  Ori  do  not  show  any  sign  of 
silicate  grains,  similar  observations  of  targets  that  belong  to  the 
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intermediate  class  group  of  Lorenz-Martins  &  Pompeia  (exhibit¬ 
ing  APO3  and  silicate  grains)  would  be  required  to  study  the  de¬ 
tailed  relationship  among  the  AI2O3  shell,  the  silicate  shell,  and 
the  photospheric  pulsation. 

Both  our  MIDI  data  and  the  images  of  the  SiO  maser  shell 
do  not  show  any  signature  of  a  globally  asymmetric  distribution 
of  dust  or  gas  in  the  case  of  S  Ori  within  our  uncertainties.  In 
general,  asymmetric  distributions  of  gas  and  dust  are  to  be  ex¬ 
pected  for  some  Mira  stars  (e.g.,  Weigelt  et  al.  1996;  Boboltz 
&  Diamond  2005;  Tatebe  et  al.  2006).  A  detection  of  asymmet¬ 
ric  structures  can  be  performed  at  the  VLTI  by  sequentially  us¬ 
ing  baselines  of  similar  length  but  different  orientation  close  in 
time/stellar  phase  to  MIDI  and  by  using  the  closure-phase  in¬ 
strument  AMBER  with  baseline  triangles  comprised  of  baselines 
with  different  orientation. 

The  relative  positions  of  star  and  SiO  maser  shells  in  our 
current  study  is  unknown.  Our  comparison  of  the  extension  of 
the  SiO  maser  shells  to  the  continuum  photospheric  diameter  is 
based  on  the  assumption  that  the  center  of  the  maser  spot  dis¬ 
tribution  coincides  with  the  center  of  the  star.  This  assumption 
is  the  most  likely  scenario  for  stars  that  do  not  show  any  sign 
of  asymmetric  dust  and  gas  distributions,  such  as  S  Ori.  In  prin¬ 
ciple,  the  relative  position  of  the  maser  shells  could  be  related 
relative  to  the  radio  photosphere,  which  was  detected  by  Reid 
&  Menten  (1997)  for  a  sample  of  6  Mira  and  semi-regular  vari¬ 
ables.  Additional  measurements  of  PRO  and  OH  masers  relative 
to  the  SiO  masers  could  give  additional  information  on  the  struc¬ 
ture  and  kinematics  of  the  more  extended  CSE. 
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